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Abstract 

We study in detail the mechanism of baryon and lepton asymmetry generation in the frame- 
work of the i/MSM (an extension of the Standard Model by three singlet fermions with masses 
smaller than the electroweak scale) . We elucidate the issue of CP-violation in the model and 
define the phase relevant for baryogenesis. We clarify the question of quantum-mechanical 
coherence, essential for the lepton asymmetry generation in singlet fermion oscillations and 
compute the relevant damping rates. The range of masses and couplings of singlet leptons 
which can lead to successful baryogenesis is determined. The conditions which ensure sur- 
vival of primordial (existing above the electroweak temperatures) asymmetries in different 
leptonic numbers are analysed. We address the question whether CP-violating reactions 
with lepton number non-conservation can produce leptonic asymmetry below the sphaleron 
freeze-out temperature. This asymmetry, if created, leads to resonant production of dark 
matter sterile neutrinos. We show that the requirement that a significant lepton asymmetry 
be produced puts stringent constraints on the properties of a pair of nearly degenerate singlet 
fermions, which can be tested in accelerator experiments. In this region of parameters the 
z/MSM provides a common mechanism for production of baryonic matter and dark matter in 
the universe. We analyse different fine-tunings of the model and discuss possible symmetries 
of the i/MSM Lagrangian that can lead to them. 



1. Introduction 



This paper is a continuation of the works [T] - [13] addressing the cosmological and phe- 
nomenological consequences of the i/MSM. The i/MSM (neutrino Minimal Standard Model) 
is a renormalizable extension of the Standard Model (SM) by three light singlet fermions - 
right-handed, or sterile neutrinos. Amazingly enough, this simple theory allows to solve in a 
unified way four observational problems of the SM [21 US] . It leads to neutrino masses and 
thus gives rise to neutrino oscillations, absent in the SM. It provides a candidate for dark 
matter particle in the form of a long-lived sterile neutrincQ, discussed already in \TQ\ [T7] . 
It allows for baryon asymmetry generation due to coherent oscillation of the other two singlet 
fermions [21 [18] and electroweak anomalous fermion number non-conservation at high tem- 
peratures [19] , associated with sphalerons [20] . A non-minimal coupling of the Higgs field to 
gravity would lead to inflation consistent with cosmological observations |13j . 

Let us recall the essential features of the z^MSM. The lightest sterile neutrino A^i plays 
the role of the dark matter particle. It should have a mass above 0.3 keV [211 [22| [23] (the 
most conservative Tremaine-Gunn bound), following from observations of rotational curves 
of dwarf galaxies. It is practically decoupled from other fields of the Standard Model (its 
Yukawa coupling to the Higgs and active neutrino must be smaller than 10~^^, as follows from 
the requirement that the mass density of the sterile neutrinos produced in the early universe 
does not exceed the dark matter abundance [H [9] and from astrophysical X-ray constraints 
[3]). Because of the very weak coupling, it does not contribute to the mass matrix for active 
neutrinos [UllK?] and does not play a role in baryogenesis. The two other singlet fermions, N2 
and N3, have masses above 140 MeV (the constraint is coming from accelerator experiments 
combined with BBN bounds [TT], see also [2^). These particles must be nearly degenerate 
in mass, to ensure coherent CP-violating oscillations leading to baryon asymmetry of the 
universe [21 118] . 

In this paper we are going to study several related issues. The first one is an elaboration of 
the mechanism of baryogenesis in singlet fermion oscillations. Though the master equations 
for leptogenesis in the i/MSM have already been written in [2], they were only analysed in 
the symmetric phase of the electroweak theory and in the regime where all singlet fermions 
are out of thermal equilibrium. We would like to answer the following questions. What is the 
CP-violating phase driving the baryogenesis? What happens with lepton asymmetry when 
these particles thermalize? When is the coherence of their quantum- mechanical oscillations, 
crucial for the resonant lepton asymmetry production lost? Are the effects of electroweak 
symmetry breaking essential for leptogenesis? What are the rates of non-conservation of 
different leptonic flavours in the z^MSM? What is the time evolution of deviations from 
thermal equilibrium due to singlet fermions, essential for leptogenesis? Can the primordial 
lepton asymmetries be protected from erasure in the i/MSM? 

The second problem is a phenomenological one. The parameter space of the i^MSM which 

number of interesting astrophysical applications of keV scale sterile neutrinos can be found in [14j . 
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can lead to the observed baryon asymmetry has never been explored in detail. Therefore, 
we would like to improve the existing cosmological constraints on masses and couplings of 
singlet fermions, which could be helpful for their experimental search. 

The third question we address in this paper is: Can we have large lepton asymmetries well 
below the electroweak scale? The motivation for this consideration is the following. In [9] 
we computed the abundance of dark matter sterile neutrino in the Dodelson-Widrow (DW) 
scenario [l5] (for earlier works see [171 [251 ESI EZ] ) . This scenario assumes that 

(i) no sterile neutrinos existed at temperatures above 1 GeV; 

(ii) the only interactions the sterile neutrinos have are those with the ordinary neutrinos; 

(iii) the universe was (leptonic) charge symmetric at temperatures below 1 GeV. 

The result was compared with two types of astrophysical bounds. The first one deals with 
X-ray observations of diffuse X-ray background of our and distant galaxies and Milky- Way 
satellites |17j . |28j - |41j and gives an upper limit on the mixing angle of dark matter sterile 
neutrino as a function of its mass. The second bound limits the free streaming length of the 
dark matter particle from observation of Lyman-a clouds [42] - [46j . The prediction, even 
with the largest uncertainties resulting from poor knowledge of QCD dynamics at the epoch 
of the quark-hadron crossover, is in conflict with astrophysical bounds. This rules out the 
DW mechanism as a source of sterile neutrino production, if one takes for granted that the 
results of [HI [l5] are robust. If the weaker, but more conservative Tremaine-Gunn bound is 
applied, then the DW mechanism can account for sterile neutrino dark matter in the universe, 
provided the mass of sterile neutrino is below 3.5 keV (the most conservative bound is 6 keV, 
see |9J). 

Since there are three essential assumptions involved, these considerations force to challenge 
one or more of them. As was found in [6j, (i) and (ii) are not valid in an extension of the 
i/MSM by a light scalar singlelEI, interacting with the dark matter sterile neutrino. The decays 
of this scalar field provide an efficient mechanism for the production of dark matter particles, 
which is in perfect agreement with all astrophysical constraints. 

In [16] it was shown that the assumption (iii) is also crucial. Namely, Shi and Fuller (SF) 
demonstrated that large lepton asymmetries can boost the transitions between active and 
sterile neutrinos leading to a possibility of resonant creation of dark matter sterile neutrinos, 
satisfying both the Lyman-a and X-ray constraints even if (i) and (ii) are correct. Qualita- 
tively, the presence of a lepton asymmetry changes the dispersion relation for active neutrinos 
in a way that it intersects with the dispersion relation for the sterile neutrino at some par- 
ticular momentum. The level crossing leads to a transfer of the leptonic excess in active 
neutrinos to the sterile ones, so that the dark matter abundance is roughly proportional to 
the lepton asymmetry. 

However, for the mechanism to work, the required lepton asymmetry must exist at tem- 
peratures 0(1) GeV and must be much larger than the baryon asymmetry, ^^^3 x 10^, 

^In ref. [6] this scalar boson was playing the role of the inflaton. The same mechanism was used in a similar 
model for another choice of parameters in [471 148] . 
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[16\ I49j . At the same time, in the majority of the models with baryon and lepton number 
violation, proposed so far, the lepton asymmetry is of the same order of magnitude as the 
baryon asymmetry. The reason is that the source of lepton number violation is associated 
with an energy scale which is of the order or greater than the electroweak scale Myy. For 
example, in Grand Unified Theories the baryon and lepton numbers are broken at the scale 
of the order of Mqut ~ 10^^ GeV, in see-saw models the masses of singlet Majorana fermions 
are, as a rule, greater than 10^ GeV. In the SM (or in its supersymmetric extension), the 
breaking of lepton and baryon numbers is related to anomaly, without any other violation 
terms. In these models baryogenesis takes place at temperatures T > Myy, and the equilib- 
rium character of sphaleron processes ensures the relation B = aL |50) - [54) . where o" is a 
coefficient of the order of one, depending on the particle content of the Standard Model or 
its extension^. A lepton asymmetry so small is irrelevant for the Big Bang Nucleosynthesis 
(BBN) and for dark matter sterile neutrino production. In conclusion, the existence of large 
lepton asymmetries seems to be very unlikely, if not impossible. 

In this paper we will show that this is not necessarily the case in the i^MSM. Indeed, the 
z/MSM is very different from the models mentioned above. In particular, the energy scale 
of the breaking of lepton number L, existing due to Majorana neutrino masses of singlet 
fermions, is small (below the electroweak scale), whereas the only source for baryon number 
(B) violation is the electroweak chiral anomaly. We will see that these facts change the 
situation so that the generation of (large) leptonic asymmetries becomes possible. Basically, 
the baryon asymmetry of the universe is related to the lepton asymmetry at the temperature 
of the sphaleron freeze-out Tew, and the lepton asymmetry generation below Tew leaves 
no trace on baryon asymmetry. The requirement that large enough lepton asymmetry is 
generated below the electroweak scale puts a number of stringent constraints on the properties 
of the singlet fermions, which can be tested in a number of accelerator experiments, discussed 

mm- 

Motivated by the fact that large low temperature lepton asymmetries can be a consequence 
of the z^MSM, in an accompanying paper [19] we reanalyze the SF mechanism for production 
of dark matter sterile neutrinos in a charge asymmetric medium. This can be rigorously done 
with the use of the formalism of [8] that allows the computation of the abundance of dark 
matter neutrinos from first principles of statistical mechanics and quantum field theory. In 
particular, we find the spectra of dark matter neutrinos which can be used in warm dark 
matter simulations, in the subsequent Lyman-a analysis and for the study of core profiles of 
dwarf spheroidal galaxies. In [49j we also establish a lower bound on the leptonic asymmetry 
A = AL/L = {riL — nL)/{nL + nLj > 2x 10~^ which is needed to make the SF mechanism 
for sterile neutrino production consistent with X-ray and Lyman-a observations (here n l and 
fii are the number densities of leptons and anti-leptons correspondingly; other conventions 

^ A breakdown of the relation L ~ _B may happen in Affleck-Dine baryogenesis [55] , if it takes place below the 
electroweak scale and if because of some reason the decay of squark-slepton condensate produces considerably 
more leptons than quarks. 
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to characterize the presence of a non-zero lepton asymmetry are described in Appendix A of 
|49j ) ■ Only this result from [l9] will be used in the present paper. 

Our findings concerning the parameter-space of the i^MSM, leading to correct baryon asym- 
metry and dark matter abundance, leads us to the fourth question we address in this paper. 
Namely, we will identify different fine-tunings between the parameters of the model, required 
for its phenomenological success, and discuss their possible origins. 

Throughout the paper we will assume the validity of the standard Big Bang theory below 
temperatures of order 1 TeV and that the only relevant degrees of freedom are those of the 
z^MSM (i.e. of the Standard Model plus three singlet fermions). After all, one of the strong 
motivations for considering the i^MSM as a theory providing the physics beyond the SM is 
the possibility to explain neutrino oscillations, dark matter, inflation and baryon asymmetry 
in the framework of a minimal model, and the creation of a baryon asymmetry requires the 
presence of temperatures above the electroweak scale. We also assume that at temperatures 
well above the electroweak scale the concentrations of all singlet fermions were zero and 
that the universe was lepton and baryon charge symmetric at this time. This type of initial 
conditions may arise in the j^MSM where inflaton is associated with the SM Higgs field [13j. 
Note also that a number of calculations in this work are on the level of approximate estimates 
which are valid within a factor of a few and thus may be refined. However, since even this 
analysis happened to be rather involved, we prefer to postpone the detailed study until the 
z/MSM gains some direct experimental support. 

The paper is organised as follows. In Sec. [2] we review the Lagrangian of the z^MSM and fix 
the notation. We also discuss different contributions to the mass difference of singlet fermions, 
essential for baryogenesis and formulate several possible scenarios for its value. In Sec. [3] we 
analyse the structure of CP-violation in the model and identify the the CP-violating phase 
that drives baryogenesis. In Sec. [4] we set up the master equations for analysis of kinetics 
of leptogenesis. In Sec. [5] we analyse CP-odd deviations from thermal equilibrium and in 
Sec. [6] CP-even perturbations. In Sec. [7] we examine in the mechanism of leptogenesis via 
singlet fermion oscillations. We derive constraints on the masses and couplings of neutral 
leptons from the requirement that the produced baryon asymmetry has the observed value 
and analyse the question whether large lepton asymmetries, which can boost the dark matter 
production, can be generated below the electroweak scale. We also determine the parameters 
of the model which allow for the survival of primordial lepton asymmetries to low temper- 
atures and are consistent with the observed baryon asymmetry. In Sec. [8] we discuss the 
fine-tunings and possible symmetries of the z/MSM and speculate on the origin of the i/MSM 
Lagrangian. Sec. [9] is conclusions, where we summarize the results. 
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2. The z/MSM and constraints on its parameters 



For our aim it is convenient to use the Lagrangian of the i^MSM in the parametrisation of 
Ref. [7]: 

J^uMSM = Cq+AC, (2.1) 

Co = CsM + Niid^rNi-{K2LaN2^ + MN2''Ns + h.c.), 

AC = -K^L^N^^-KiL^Ni^-^^ NfNj + \i.c., 

where Nj are the right-handed singlet leptons, $ and (a = e, /i, r) are the Higgs and 
lepton doublets respectively, /i is a matrix of Yukawa coupling constants, M is the common 
mass of two heavy neutral fermions, different elements of AM/j ^ M provide a mass to the 
lightest sterile neutrino A^i, responsible for dark matter (Mi ~ AMn) and produce the small 
splitting of the masses of N2 and A^3, AM23 = 0, $i = ^ij^]-, and M is taken to be real. 

Another parametrization of the same Lagrangian is related to the mass basis of Majorana 
neutrinos. We write 

M = UXU^ , (2.2) 

where Aiij = M6j25j3 + M5i^5j2 + AM/j, A is a diagonal matrix with positive values A/, 
and [/ is a unitary matrix. The values A^ are nothing but eigenvalues of the hermitean matrix 
M.M.'^ which can be diagonalized with the help of U , Ai ~ Mi. 

Yet another possibility is to usG the bOiSis in which the matrix of Yukciwa couphngs h(y^i is 
diagonal, 

h = KLfdKl (2.3) 

with fd = diag(/i, /2, /s). Definitions of the matrices Kl and Kr can be found in [2]. 

The Yukawa coupling constants of the dark matter neutrino A'^i are strongly bounded by 
cosmological considerations [1] and by X-ray observations [HEI]: J2a I l^aiUnl'^ < 10"^'^. 
As has been demonstrated in [1], the contribution of the dark matter sterile neutrino to the 
masses of active neutrinos via the see-saw formula 

m.u = -^^^^^^ (2.4) 

is much smaller than the solar neutrino mass difference, and can thus safely be neglected 
(here M^j^i = haiv and v = 174 GeV is the vacuum expectation value of the Higgs field). 
Therefore, we set h^i = in the following and omit A^i from the Lagrangian. Note that eq. 
(j2.4p implies that the mass of one of the active neutrinos is much smaller [Ij than the solar 
mass difference ~ 0.01 eV and can be put to zero in what follows. 

The values of the Yukawa coupling constants ha2 and h^s are further constrained by the 
requirement that the i^MSM must describe the observed pattern of neutrino masses and 
mixings. The following relation must hold: 

r ,.2 



M 



N 



(2.5) 
IJ 
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where M,y is the mass matrix of active neutrinos, and we denoted hy A4n the 2x2 mass matrix 
of the second and third singlet fermion, [A^atJ/j = Aijj for I,J = 2,3. This formula can be 
simplified further by noting that -^3 rnust be highly degenerate in mass in order to ensure 
successful baryogenesis [2l|T8]. In fact, any non-zero mass difference that still remains (it will 
be discussed later) is inessential for discussion of masses and mixings of active neutrinos and 
can be ignored [7j. We have, therefore, 

[Mu]af3 = -J^ {ha2hp3 + /ia3/i/32) • (2.6) 

As was shown in ref. 0, this simplified situation allows to determine Yukawa coupling con- 
stants from the mass matrix of active neutrinos up to rescaling ha2 ha2/z, ha3 zhas, 
where z is an arbitrary complex number. In addition, one can solve for the active neutrino 
masses explicitly: 

m G {0,t;2[F2F3 ± \h^h\23]/M} , (2.7) 

where = [h^hla. This leads to two qualitatively different cases, namely the "normal 
hierarchy", mi = 0, m2 = msoi, ^3 = Jn-atrnj and the "inverted hierarchy", mi ~ m2 ~ 
m-atm, "I3 = 0. Here mgoi = y^Am^^^, m^tm = y^Am^, and Am^„^ ~ 8.0 x 10"^ eV^, 
Am^^^ — 2.5 X 10~^ eV^ [56]. Normal hierarchy corresponds to the case |/i^/i|23 ~ F2F3, and 
the inverted hierarchy to the case |/i^/i|23 ^ F2F3. From here it follows that 

2F2F3V^/M ~ Kmatm , (2.8) 

where k = 1 (2) for normal (inverted) hierarchy, matm ~ 0.05 eV. If F3 is taken to be very 
small, F = F2 is required to be large to keep the atmospheric mass difference in the right 
place. The ratio of Yukawa couplings F2 and F3 will play an important role in what follows 
and is denoted by e, 

^ = 5 • (2.9) 



Then 

^2 _ KMrUgtm 

2ev 



i72 = ^^^^^^1^. (2.10) 



In the limit e ^ 0, AM/j the Lagrangian acquires the global leptonic U(l) symmetry 
[7] which guarantees the degeneracy of the A'^2 and N3 of singlet fermions (necessary for 
baryogenesis) , absence of mass for A'^i and absence of interactions of A^i with active fermions 
(providing thus an approximate description of the required parameters of the dark matter 
sterile neutrino). This symmetry can be made explicit by introducing the 4-component Dirac 
spinor ^ = N2 + iVJ unifying a pair of two degenerate Majorana fermions. 

For the discussion of the baryon and lepton asymmetries of the universe an essential pa- 
rameter is the mass difference 6M between the mass eigenstates of the two heaviest neutrinos 
[21[T2]. Indeed, successful baryogenesis can take place provided the mass difference is small 
enough. In the theory defined by the Lagrangian ()2.1I) there are two sources for the mass 
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difference: the first one is related to the Majorana mass matrix, and the second one is due to 
the Higgs vacuum expectation value and Yukawa couplings to active fermionic flavours. The 
mass square difference of the physical states at the leading order of perturbation theory with 
respect to Yukawa couplings and AMjj is given by 

I 2 1 

»/ = L_l, (2.11) 

where 

= 2{h^h)23V^ + M{AM^2 + AM33) . (2.12) 

The one- loop corrections to this result are of the order of ^I^K . 

One can distinguish three essentially different situations, depending on the relative impor- 
tance of the mass difference induced by the Higgs field and the difference associated with 
Majorana masses. 

In the first one AA = A3 — A2 is negligibly small so that the mass difference is entirely 
due to the Higgs condensate. One can easily find in this case from ()2.7l2.12p that the mass 
difference of heavy neutrinos is the same as that of active neutrinos, 

6M = Aniy , (2.13) 

and is given by 6M = rriatm ~ f^soi — 0.04 eV for the case of normal hierarchy or 6M ~ 
Am2^i/2matm ^ 8 x 10"'' eV for the case of inverted. Quite amazingly, these mass differences 
are roughly those at which the production of baryon asymmetry is extremized, M5M ~ 
where Mpi is the Planck mass, Mpi = 1.22 x 10^^ GeV (see below for a more detailed 
discussion). We will refer to this situation as Scenario I for singlet fermion mass difference. 

The second option is when the mass differences coming from two different sources are of 
the same order of magnitude. We will call this choice of parameters Scenario II. An extreme 
case is an exact compensation of the two leading contributions, which would allow to have a 
(low temperature) mass difference be much smaller than the active neutrino mass difference. 
To realize this fine tuning, the following condition is required to hold: 

5M < Amy . (2.14) 

Though this possibility may be considered to be bizarre (as it requires that contributions of 
seemingly different nature are exactly the same in magnitude and different in sign) it must 
not be discarded since the origin of different terms in the z^MSM is unknown, so that eq. 
(j2.14p could be a consequence of some underlying (Planck scale?) physics. We will call this 
option Scenario Ila and discuss this fine-tuning in more detail in Sec. [H 

Yet another possibility is when AA ^ Arrii, so that the mass difference is entirely due to 
the Majorana masses. We will refer to this possibility as Scenario III. 

For all three cases the mass eigenstates are related to the fields by a rotation with the 
maximal angle 7r/4 (up to complex phases and corrections ~ 0{M£)/M)). 
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3. The structure of CP violation in the z/MSM 



The generic Lagrangian (j2.ip contains a number of new physical parameters (18) in com- 
parison with the Standard Model. They can be counted as follows: 3 Majorana masses of 
singlet fermions, 3 Dirac masses, 6 mixing angles and 6 CP-violating phases. As we have 
already mentioned, one of the singlet fermions, playing the role of dark matter, has very 
small Yukawa couplings, and thus is irrelevant for baryogenesis and for active neutrino mix- 
ing matrix. Thus, the number of parameters of the j^MSM, responsible for physics of heavier 
singlet leptons is smaller. The aim of this Section is to identify these parameters and to find 
the CP-violating phase, relevant for baryogenesis. 

After dropping the dark matter sterile neutrino A^i from the Lagrangian (j2.ip the theory 
contains 11 new parameters in comparison with the SM. These are 2 Majorana masses, 2 
Dirac masses, 4 mixing angles and 3 CP-violating phases. In principle, all these 11 physical 
parameters of the I'MSM can be determined experimentally by the detailed study of decays 
of the singlet fermions. At the same time, the mass matrix of active neutrinos in this case 
(note that one of neutrinos is massless in this approximation) depends on 7 parameter^: 3 
mixing angles 9i2, 623 and ^13, one Dirac phase (p, one Majorana phase (a in the normal 
hierarchy case and the combination C, = {a — 0)/2 in the inverted hierarchy) and two masses 
(m2, for normal hierarchy, and mi, m2 for inverted). So, the 11 new parameters of the 
z^MSM can be mapped to 7 parameters of the active neutrino mass matrix plus 4 extra ones. 
It is convenient to select these 4 parameters as follows. The first one is the average Majorana 
mass of singlet fermions A^2 and N^, 

M = ^^±^. (3.1) 

The second is related to the diagonal elements of the Majorana mass matrix of singlet 
fermions, which after phase transformations without loss of generality can be written as 

AMm = ^^y^ = AM22 = AM33 . (3.2) 

The third one is the parameter e defined in ()2.9p . and the fourth is an extra CP-violating 
phase ?7, associated with it (see below). 

With these notations, the relevant part of the i^MSM Lagrangian is: 



L 



singlet 



1 

KMniatm \ 2 



2v 



2 



^L2iV2 + Ve^LsNs 



MN2^N3 - ^AMMiN2^N2 + iV3"iV3) + h.c. , (3.3) 



where L2 and L3 are the combinations of Le, Lu and Lj- 



-^2 = ^ , -t-S = ^ • (3.4) 

^3 



*We use the notations and parametrisation of ref. [56] 
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We stress that the relation between L2,3 and -^e,^,T is not unitary, in general. In fact, there 
are 4 different relations between L2,3 and leptonic flavours for each type of hierarchy [7J (some 
of them may be equivalent to each other after phase redefinition of leptonic flavours), leading 
to one and the same active neutrino mixing matrix. We will present and analyse below only 
one of them for each hierarchy, others can be treated in a similar way. 

To simplify the analysis, we consider the case (suggested by experiments) when the angle 
6i3 and deviation of 6*23 from its maximal value, 8623 = O23 
the notations 



J are small. Let us introduce 



Di = 6923 cos 012 + 013 sin 0i2e^<^ , D2 = 8623 cos 9i2 - ^13 sin e^e'''' 
D3 = (56'23 sin 6*12 + % cos 6*126*'^ , L»4 = 56*23 sin 6^12 - 6*13 cos 6li2e*'^ 

Then, for the normal hierarchy case we can write: 

L„ - Lt 



L2 
L3 



+ai(l + zi + Z2) 



-ai(l - zi- 
where the main terms are 

ai 

03 
tan p 

The corrections are given by 

zi 

Z2 



Z2) 



Ln + Lt 

--j= 1 a2(l + Z2- Zi)Le + 03 ^ 

L„-Lt , at , L^ + Lt 

02(1 - Z2 + Zi)Le + a3-^ 



V2 



V2 



^g-i(a+<^) sin p COS 012 



le smpsmt/12 
cosp , 

Im2 




0.4 . 



cot p 



' sin 2012 ' 

^•n*ei(-+</')^^ - e-^('^+<^)L»itanp 
^ sin 2012 



For the inverted hierarchy the corresponding equations are: 

L 



L2 
L3 



+ie-'ni{l + t\-tl)- 



V2 



Lt , , , ^ , .Ln + Lt 

+ b2{l + ti+t2)Le + {z3-t3)- ^ 



V2 



+ bUl -h+ t2)L, + {Z3 + t3)-^^ + 



V2 



V2 



where the main terms are 



bi 
62 



1 
1 

72 



cos 0126 *^ + Z sin 0126"''*'' 



cos 0126+*^ + i sin 0i2e 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



(3.10) 
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and the corrections are given by 

2i cos 2^ sin 2^12 



ti = 5„ 



t2 = 6i 



3 + cos 4C + 2 sin^ 2C cos 46*12 
1 1 



.2 1 + e-4< tan2 6*12 



(3.11) 



where 

'm2 — mi \ Arn^ 



fJ- ~ 8 X 10"^ . (3.12) 



m2 + miJ 4.Aml^^ 

In general, baryon asymmetry of the universe in the z^MSM may depend on all three 
CP-violating phases described above. However, in a specific limit, when all charged lepton 
Yukawa couplings are the same, the baryogenesis is driven by a single phase, which we identify 
below. This can be done on general grounds and does not require complicated computation. 

If the Yukawas in the charged sector are the same, one can choose a basis for leptonic 
doublets, in which interactions of singlet fermions have a simple form 

{f2l2N2 + fskNs + f23l2N3)^ , (3.13) 

where ^1,2,3 are related to Le^^^r by some unitary transformation. Now, by the phase redefi- 
nition of I2 and /s the constants /2 and /s can be made real. Then, in this parametrization, 
the CP-violation effects must be proportional to the complex phase of the coupling /23, they 
must vanish in the limit /23 — > 0. It is easy to see that 

A3 = . (3.14) 

Now, with the use of eqns. (j3.7l3.8p we get for the normal hierarchy : 

/23 = eFe^^ cos 2p [1 — 2i tan p ((56*23 cos(a + (/>) cos O12 + ^13 cos a sin ^12)] . (3.15) 
The similar expression for the inverted hierarchy is obtained with the use of (|3.10|3.1ip : 



/23 = -^Fe'"^ 



Sinv + \it*3- 4) (*3 + Zs] 



(3.16) 



Yet another parameter in (j3.3p which is important for the issue of CP-violation is the 
mass splitting AMj\,/. Indeed, in the limit AMm = (the extreme case of Scenario I), 
the Lagrangian (j3.13p acquires the global U(l) symmetry, so that the phases of A'^2,3 cannot 
be fixed anymore by the mass terms, and the CP phase of /23 in the interaction ()3.13p can 
be rotated away. This means that the measure of CP-violation, relevant for baryogenesis 
for the case when all charged lepton Yukawa couplings are the same, can be conveniently 
parametrised as 

6cp = e sin(arg 723) sin 9 , (3.17) 
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where 

tane = ^. (3.18) 

In reality the charged lepton Yukawa coupUngs are different, and the structure of CP- 
breaking relevant for baryogenesis is more involved. It has been found in [2] for the case 
when all reactions of singlet fermions are out of thermal equilibrium, see eq. (29) of that 
papei^l. For convenience, we present it here in notations of our work, factoring out the largest 
Yukawa coupling F: 



(3.19) 

With the use of relations (|3.4|3.6|3.9p . the expression (j3.19p can be rewritten through the 
parameters of the neutrino mixing matrix, the value of e and phase rj. The corresponding 
relations are not very illuminating, so that we just summarize the main qualitative features 
of (|3T9D . 

(i) The sign of baryon asymmetry cannot be found even if the active neutrino mixing matrix 
is completely known. 

(ii) For small e, 6cp oc e, similar to eq. (|3.17p . 

(iii) Baryon asymmetry is non-zero even if ^13 = and 8623 = 0. In other words, the details 
of the active neutrino mass matrix have little influence on baryogenesis, if these angles are 
small. 

(iv) Baryon asymmetry does not vanish even if tan 9 = 0, i.e. in the Scenario I for singlet 
fermion mass differences. In this case the result is determined by the parameters of the 
neutrino mixing matrix only. 

(v) In general, (Jcp 7^ for e = 1. 

(vi) For inverted hierarchy case, 6cp 7^ even if one takes a limit rrii = m2, ^13 = 0, 6623 = 0. 
The computation, leading to eq. (I3.19p . cannot be applied to the low temperature leptoge- 

nesis that occurs at temperatures in the GeV range (see Sec. [7]), because the mass of r-lepton 
is comparable with a relevant temperature and with the mass of singlet leptons A''2,3. In 
this case the charged lepton masses certainly cannot be neglected, and the structure of CP- 
violation is even richer than that above the electroweak temperature. Still, CP-violating 
effects are proportional to e, since in the limit e — > all couplings in (|3.3p can be made real. 

For numerical estimates of CP-violating effects in the paper we will assume that the relevant 
CP-violating phase is of the order of one, so that the effects are suppresses by e, writing 
explicitly this factor in the formulas. 



^The fact that this result does not depend on Yukawa coupHngs of charged fermions does not mean that 
they can be neglected. Indeed, the computation of baryon asymmetry takes into account that all reactions of 
the particles of the SM equilibrate, and the fact that charged Yukawas are non-zero and different is essential 
for kinetic description of the system with the use of eqs. (|4.3I4.4I4.5|) . 
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4. Lepton asymmetry generation: review of theoretical framework 

The detailed description of the system of singlet leptons and active fermions in the early 
universe is necessarily quite complicated. The number of relevant zero-temperature degrees 
of freedom (3 active and 3 sterile neutrinos and their antiparticles) is largqj, and the time- 
scales of different processes can vary by many orders of magnitude. Moreover, due to the 
smallness of the sterile-active Yukawa couplings the processes with singlet fermions have in 
general a coherent character, making the approach based on Boltzmann equation for particle 
concentrations useless. 

Probably, the simplest way to deal with coherent effects is to use the equation for the 
density matrix [571 ESI HSl [2]. In our case this is a 12 x 12 matrix (12 = 3x2x2 degrees of 
freedom for all active and sterile neutrino states), satisfying the kinetic equation (11) of [2]: 

i^ = [H,p]- ^^{r,p}+'-{rp,i-p} , (4.1) 

where H = p{t) + Ho + Hint is the Hermitean effective Hamiltonian incorporating the medium 
effects on neutrino propagation, p{t) is the neutrino momentum, with {p{t)) ~ 3T (we will as- 
sume that all the neutral fermion masses are much smaller than the temperature), Hq = 
(we include AM/j to Hint), T and are the Hermitean matrices associated with destruc- 
tion and production rates correspondingly, and [,]({, }) corresponds to the commutator 
(anti-commutator) 0. Following refs. |18l [2] we will use the Boltzmann statistics for esti- 
mates and replace the last term in (j4.ip by iT^. Also, following [2] we will replace by 
^{T,p'^'^} with yO^' = exp {—p/T) being an equilibrium diagonal density matrix, ensuring the 
correct approach to thermal equilibrium. After these substitutions the kinetic equation takes 
a simple form 

i^ = [H,p]-'-{r,p-p^'^} . (4.2) 

This is a relaxation time approximation for the density matrix, fairly standard one in non- 
equilibrium statistical physics. 

This equation can be simplified even further (for details see [2]) accounting for the following 
facts: 

(i) The rates of interactions between active neutrinos are much higher that the rate of the 
universe expansion. Therefore, coherent effects for active neutrinos are not essential and 
the part of the general density matrix p related to active leptonic flavours can be replaced 
by equilibrium concentrations characterised by 3 dimensionless chemical potentials pa (the 
ordinary chemical potential divided by the temperature) giving the leptonic asymmetry in 

®In fact, the number of types of particle excitations in high temperature plasma is even higher, but we will 
assume in this paper that only zero-temperature degrees of freedom are relevant. 

'^We stress that ea. (|4.ip . though it looks identical to eq. (1) of the earlier work [18], is in fact very different. 
In ref. [18] p is a 3 x 3 matrix, associated with singlet fermions only, whereas ea. (|4.H) accounts for all leptonic 
degrees of freedom. 
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each flavour. 

(ii) Active neutrinos get temperature dependent masses that are quite different from those 
of singlet fermions. Therefore, all non-diagonal elements of the density matrix involving 
simultaneously the active and sterile states can be put to zero. 

(iii) The coupling of the dark matter neutrino is so weak that it decouples from the system. 
This leaves us with the 2x2 density matrix for singlet fermions A'2 and N^, charge 

conjugated density matrix pN for corresponding antiparticles (or, to be more precise, opposite 
chirality states), and 3 chemical potentials pa- The corresponding equations can be written 
as [2]: 

= [H,pN]-'-{TN,PN-p'^ + ipian ^ (4.3) 
.d^ = [H\-p^]-'-{T*^,pN-p''}-il^an* , (4.4) 

.dpa 



dt 



iVlp^ + iTi TUpn-p''') -iT, TfipN-p''') . (4.5) 



In the equation for p^ there is no summation over a and Tg are real. The explicit expres- 
sions for the matrices describing different equilibration rates (r^r, F^, T^, T^) via Yukawa 
coupling constants can be found in [2j for the case when the temperature is higher than the 
electroweak scale. They are all related to the absorptive parts of the two point functions 
for active or sterile neutrino states and contain a square of Yukawa couplings hai- The real 
parts of the corresponding graphs together with mass squared difference between and 
determine the effective Hamiltonian H. For high temperatures T > Tew the equilibration 
processes are associated with Higgs, W and Z decays to singlet and active fermions, to cor- 
responding inverse processes, and to ti — > Ni' scattering (t is the top-quark). At smaller 
temperatures T < Tew the rates are associated with W and Z exchange and singlet-active 
mixing through the Higgs vev, see Fig. [TJ 

In the earlier work [18] the computations and qualitative discussion were based on incom- 
plete kinetic equations, which did not include the last terms in eqns. ()4.3p and (14. 4|) . as well 
as eq. (|4.5p . As was shown in [2], these terms are absolutely essential for the analysis of 
lepton asymmetry generation. Therefore, all the quantitative results for baryonic excess of 
[21 [7] and of the present paper are different from those of [18]. If some of the qualitative 
conclusions happen to be the same, we cite both papers [18] and [2] simultaneously. 

The eqs. (|4.3|4.4l4.5p . supplemented by an initial condition pn = Pn = Pa = ^ which may 
be fixed by inflation, provides a basis for the analysis of the lepton asymmetry generation [2]. 
In that paper an analytic perturbative expression for lepton asymmetry has been derived. 
The expression (33) of [2] is valid provided the following requirements are met: 
(i) All reactions, in which the singlet fermions participate, are out of thermal equilibrium 
above the sphaleron freezing temperature. In this case a straightforward perturbation theory 
on Yukawa couplings of singlet fermions can be used. If this requirement is not satisfied, 
a perturbative expansion contains the so-called secular terms, which diverge with time and 
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Figure 1: Diagrams for the processes which contribute to equilibration rates. 



require resummations. It is one of the aims of the present paper to find out what happens if 
the singlet fermions equilibrate before the sphalerons decouple. 

(ii) The mass difference between singlet fermions is sufficiently large so that the Higgs field 
contribution to it can be neglected. In other words, only Scenario III was considered. 
Naturally, we would like to extend the analysis of baryogenesis to Scenarios I, II, and Ila. 

(iii) The number of oscillations of singlet fermions, related to their mass difference (see exact 
definition in eq. (I7.14p ) is much greater than one at the time of the electroweak cross-over. In 
this case the baryogenesis occurs in the symmetric high temperature phase of the SM and the 
Higgs vacuum expectation value does not play any role. We would like to understand what 
happens if this assumption is not satisfied, which is the case, in particular, in Scenarios I, 
II, and Ila. 

To address all these questions we make a number of helpful transformations of kinetic 
equations (j4.3|4.4l4.5p . In particular, a further simplification of the system (j4.3l4.4l4.5p can 
be made under assumption that the CP-violating effects are small. Let us introduce the 
CP-odd (p-) and CP-even deviations (p+) from thermal equilibrium by writing 

p^-p-'i = 5p+ + ^^ , p^-p'^i = 6p+-^-^ . (4.6) 

and neglect in (j4.3l4.4|4.5p all terms that are of the second order in CP-odd quantities (such 
as (Pjv — T*j^)Sp-) etc. In this approximation one can decouple the equations for the CP-even 
deviations p+ and get 

= [Re H,5p+] - '-{Re Tn,Sp+} , (4.7) 
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with the initial condition 
mation have the form: 



-p'^'i. The equations for the CP-odd part in this approxi- 



.d5p- 
' dt 

.djJa 



[Re H,5p-] - -^{^<s Tat, 6p-} + i/ioRe + S, 



dt 



-iTlfia + iTr ReTl5p 



+ 5a 



with zero initial conditions for 6p- and leptonic chemical potentials. Here the source terms 
5" and are proportional to CP breaking parameters and given by: 



S 



2i[Im H, 6p+] + {Im T^, Sp+} , 



-2Tr 



Im rl6p, 



(4.9) 
(4.10) 



They are only non-zero when CP-even deviations from thermal equilibrium exist, which is a 
key issue for baryogenesis and leptogenesis [59]. At the same time, if different damping rates 
in (j4.8p are all larger than the rate of the universe expansion after leptogenesis, the created 
asymmetry disappears. Therefore, to find whether baryogenesis is possible at all, one can 
study first the rates of different processes that equilibrate CP-odd and CP-even deviations 
from thermal equilibrium. This can be only skipped if all reactions with singlet fermions are 
out of thermal equilibrium, which was the case considered in [2]. If the necessary conditions 
for baryogenesis are found to be satisfied, an analysis of the CP-violating effects must follow. 

So, we will consider first the system (j4.7l4.8p neglecting all CP-violating effects. To simplify 
the notations, we will take away the symbol Re of the real part from the equations. 



5. CP-even deviations from thermal equilibrium 

5.1. High temperature singlet fermion masses and mass eigenstates 

The behaviour of the CP-even perturbations is determined by H and F^v, see eq. (|4.7p . Let 
us start from a discussion of the Hamiltonian Hint, describing the oscillations. 
The Hamiltonian Hint has the form 

Hint = 2^~~^ ' ^^'^^ 

where AM^(T) is the temperature dependent (non-diagonal) matrix of mass differences be- 
tween singlet fermions. It is determined by the zero-temperature mass difference and by real 
parts of propagator-type graphs for sterile fermions, see Fig. [5J 

There are two different temperature-dependent contributions to the mass difference. The 
bottom one is proportional to the square of the temperature dependent vacuum expectation 
value of the Higgs field v{T) whereas the top one is proportional to T"^, coming from the 
Higgs exchange. At temperatures around and below the sphaleron freezing, interesting to us. 
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Figure 2: "Soft" contribution to the mass difference of singlet fermions coming from electroweak 
spontaneous symmetry breaking (lower panel) and from radiative correction (upper panel). Non-zero 
temperature neutrino propagator has to be used. 



the contribution related to the Higgs vev dominates because of the usual loop suppression 
and since v{T) > T. We get for the high temperature case <C T^: 



AM\T)jj 











v^{T)haih*^j- 



2bp 



(5.2) 



m*\T) J ^ M2 + 26p 

where m?{T) is determined by eq. (I2.12|) with the replacement v v{T). The function h is 
defined by the active neutrino propagator l/(j/ + ^q): ^ = + 6^/, where u is 4- vector of 
the medium. The function b in different limits is given by |601 161|: 



(2 + COS^ 6*14/) 



TTCHV \ ' " I 360 

This is the so-called potential contribution to active neutrino dispersion in the medium. 
From (j5.2p we get for the temperature-dependent mass difference 6M{T) 



6M{T) 



2M 



(/it/i) 



23 



m?{T) 



(5.4) 



For small e and high temperatures this gives in Scenarios I, II: 

5M atm_V/ ^ 

The temperature-dependent contribution to the mass difference is suppressed in comparison 
with the zero-temperature one at 2bp <^ A'P. Taking for an estimate the typical momentum 
of a particle in high temperature plasma p ~ 3T one finds that this inequality is satisfied at 



T < X, 



pot 



f M \ 



1/3 



GeV 



(5.6) 
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and that for these temperatures and e ^ 1 the mass difference is 



SMiT) ^ ^ {^^y + SMHO) . (5.7) 

Note that depending on parameters 5M{T) can go through zero at some particular tem- 
perature, leading to level crossing and to the resonant production of lepton asymmetry, see 
Sec. [71 As follows from eq. (|5.4p . this may only happen at T > Tpot if e = 1. 

Let us discuss the high temperature mass eigenstates for three different scenarios of the 
singlet fermion mass differences (see the end of Sec. [2] for definition). 

Scenario I. In this case AA = and the mass difference comes entirely from the interaction 
with the Higgs field. 

For e ^ 1 the high temperature mass eigenstates N'2 and A''3" are close to and N^, 

~ cos/3oiV2 + sin/5oA^3 , 

~ cos/3oiV3-sin/3oA^2 , (5.8) 

where 

whereas for e ~ 1 they represent the mixing of and N^, with the angle of the order of 1. 
With the use of eq. (j2.8p the ratio of Yukawa couplings that appear in (j5.9p can be written 
as 

eArriy | 1 Normal hierarchy , , 

Po — — e < g ^ (5.10) 

i^n^atm I 8 X 10 Inverted hierarchy . 

It is not difficult to see that for temperatures below Tpot the mixing angle (3 gets modified, 

/3o-/3^/3o(^l + ^j , (5.11) 

leading to ^M) for bp > (this expression is valid provided /3 < 1). Therefore, the 
temperature Tp at which the mixing (3 is of the order of one is given by 

r,.4T,„,.16Gev(l^«)*(JL)'. ,5,12) 



This is derived with the use of eq. (j5.3p for T <^ Mw and e ^ 1. 

Scenario II. In this case both terms in eq. ()5.2p have the same order of magnitude and 
the mixing angle (3 is in general of the order of one. It goes to the zero-temperature value 
7r/4 at temperatures below T^, see eq. (|5.12p . For the Scenario Ila at T > the mass 
difference is of the order of matm/^ (see eq. (15. 5p ) and is much smaller than Am,y at lower 
temperatures, see eq. (|5.7p . Note that at T > Tew ^{T) / v and, therefore, m{T) / even 
if m(0) = 0. 

Scenario III. In this case the mass difference comes entirely from the tree Majorana mass, 
and the high temperature mixing angle is always close to 7r/4. 
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5.2. Damping of CP-even perturbations 



Let us turn now to the part of equation (14. 7j) describing creation and destruction of singlet 
fermions. To find Ftv we note that for initial condition p = 0, ^ = one gets for sufficiently 
small times from eq. (|4.7p that 

dp 



(5.13) 



At the same time, for the time scales smaller than 1/Tn but much larger than microscopic 



time scales such as l/T, the derivative ^ can be found from first principles of statistical 
mechanics and quantum field theory as described in [HI [9]. Therefore, we can use the methods 
of these papers to define the rate Ftv- 

First we note that for the computation of F the term ^Mii. Nj'^Nj can be neglected (this 
term, however, must be kept in Hint, as has been done above). In this case the fields A''2 and 
A^3 can be unified in one Dirac spinor as ^' = + -^l- As usual, ^ can be decomposed in 
creation and annihilation operators as 

d^i 



:U[p,s)e 



-iP-x 



+bi 



AP-x 



where the spinors u, v satisfy the completeness relations 

^ ^(p, s)n(p, s)=i) +M, ^ t;(p, s)v(p, s) 



i> -M 



(5.14) 



(5.15) 



The operators dq^s and 6q^s are the creation operators of a singlet fermion and an anti- 
fermion A^3 with momentum q, and helicity state s. These operators are normalized as 

{«p,.,<,J = 'J^'np-q)5st, (5.16) 

and V is the volume of the system. The density matrix pAf is associated with operators 

1 / dq, + dq ^ «i 



Now, repeating literally the discussion of the Section 2 of ref. [8] we arrive at 



ri/(q) 



(5.17) 



(5.18) 



where p is the spectral function defined in Appendix B of ref. [8j, and matrices H/j are given 

by 



v''{T) [\K2?Pu{v) + \K-i?P^{-v) 

h*^^ha3v\T)Pl{p) , 

v\T) \\K2\^P,{p) + \ha3\^PS{-p) 



(5.19) 
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The spin operators Pu = u(p, +)n(p, +) and Py = v{p, +)t;(p, +) are 

Pu = l^Pul\ Pi = ^^Py^\ u' = v, v' = u, p=|p|. (5.20) 

Let us discuss the structure of the matrix Fj^ in more detail. We will be interested in a 
total (integrated over momenta) rate appearing in eqs. (j4.3l4.4l4.5p : 

Tn{T,M) = ^J d^qTi,'{q) . (5.21) 



Then the structure of F/v is: 



, ^_ R{T,M) + e'RMiT,M) ih^R(T,M) 
^ Fo\ i^^R{T,M) e^R{T,M) + RM{T,M) ' ' 

where i<o = 2 x 10~^ is a convenient normalisation constant, and 

i?(T, M) = J d^qT%\^)U=o (5.23) 

can be called the rate of the singlet fermion production at F = F^. The quantity 

^0 f jsQ - Qrii, 



Rm{T, M) = ^ J d^gl^F]v^(q)|,=o (5.24) 

vanishes in the limit M ^ and represents the rate of the processes with violation of total 
lepton number (to be defined exactly below). 

Computation of Rm{T,M) and R{T,M) is quite involved and is discussed in detail in 
Appendix A. A large number of processes, such as W, Z and Higgs decays, together with 
2^2 reactions incorporating quark and lepton initial and final states must be taken into 
account. The result of the computation is presented in Fig. [3l The vertical axis is the 
temperature T, and the horizontal axis is the temperature derivative of the yield parameter, 
defined in eq.(4.8) of [9] : 

T^ = -.{T)R{T,M), T^ = -.iT)R,dT,M), .(T) - 



dT V ^ ' V ^ -^''^v ' V / 47r2c2(r)/ie//(T)r2 ' 

(5.25) 

where Cs is a speed of sound, the temperature-time relation is given by t = Mq ~ 

Mpi/1.66^geff, and the temperature dependence of the numbers of degrees of freedom g^jf 
and heff can be taken from [9j. The combination 

1 dY , , 
T (5.26 

Yeg dT ^ ' 

is nothing but the ratio of the singlet fermion production rate to the Hubble constant. 
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Figure 3: The temperature derivative of the yield parameter related to the rate R{T, M) for the 
Higgs mass itlh ~ 200 GeV, F — Fq and different values of the singlet fermion mass (left panel). 
Right panel: the same for Rm{T, M). 



For temperatures smaller than the peak temperature the rate R(T, M) can be reasonably 
approximated by 

^i?(r,M) -SG^T^^o' , (5-27) 

where ~ 5 is a numerical constant found by fitting of the numerical result, 9q = 
is the zero-temperature mixing angle between the singlet fermion and active neutrinos. At 
the temperatures above and around the peak the suppression of the transitions due to the 
medium effects [T5l[16l[T7] becomes important. Also, the decays of the vector bosons and of 
the Higgs must be taken into account. At temperatures in the region 100 — 200 GeV the rate 
scales as R{T,M) oc l/T, while at temperatures above the peak roughly as R{T,M) oc 
In the symmetric phase of the electroweak theory, T > 250 GeV, studied previously for 
baryogenesis via singlet fermion oscillations in [21 [18], the rate scales like R{T,M) cx T. 

5.3. Time evolution of CP-even perturbations 

Having defined the mass matrix of singlet fermions and the matrix of the damping rates we 
are ready to consider the behaviour of CP-even deviations from thermal equilibrium. Let us 
choose the basis in which AM^{T)jj is diagonal: 

^E=(l' (5.28) 



£■ 



3 
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and 

Tjv = ( 1 . (5.29) 




Since for practically all temperatures AE ^ ^u>t^ '^^^ easily finds four different exponentials 
describing the time behaviour of the density matrix: 

exp[-((r22 + r33)/2±i(^2-^3))i] , exp(-r22t) and exp(-r33t) . (5.30) 

The first one corresponds to the behaviour of the off-diagonal elements of Ap_|_ and thus to 
the damping of quantum-mechanical coherence in the oscillations of singlet fermions. Two 
others represent the approach to thermal equilibrium of the diagonal elements of the density 
matrix. In general, if T22t > 1 and T^^t > 1 the system equilibrates completely. 

Let us consider again different scenarios for the singlet fermion mass matrix. 

In the Scenario I the matrices Hint and can be simultaneously diagonalised for T > 
Tpot and e <S 1 (up to the mass corrections M^/T^). Then, for e <C 1 we have two very 
different relaxation rates for the diagonal elements of the density matrix, 

r22^^i?(T,M), ^r,R{T,M) + ^RM{T,M), r, = ^ (^e^ - K^^l^l^^ , (5.31) 

whereas the rate of coherence loss is related to r22. With the use of relation (j2.7|) the 
combination of Yukawa couplings which appears in (I5.3ip can be represented as 

eMmatm I 0.36 Normal hierarchy „„n 
v'^Fq ] 1 Inverted hierarchy ' 

leading to r33/r22 oc e^. When the temperature falls down from Tpot to Tp the mixing angle 
P changes from small values ~ e to /? ~ 1, which modifies the (smaller) rate r33 as 



33 



r33 + sin2/3r22 . (5.33) 



As a result, at T < both rates are of the same order of magnitude and are related to the 
largest one r22. Of course, for e ~ 1 all damping rates have the same order of magnitude for 
all temperatures. 

In the Scenario II the matrix Hint can be diagonalised with the help of orthogonal trans- 
formation O, O^HintO = diag, characterised by the angle /3 ~ 1. In general, the rates r22 and 
r33 are of the same order. The same is also true for the Scenario III with AA > Arrii^/e, 
leading to the mixing angle /3 ~ 7r/4. For e <C 1 all the damping rates are nearly the same 
and equal to 7y^R(T, M). Qualitatively, if the rate of oscillations between strongly coupled 
singlet fermion (rate [r^*]22) and weakly interacting fermion (rate [r^*]33) is large, the ap- 
proach to thermal equilibrium is determined by the largest rate since the system spends half 
of the time in the strongly interacting state. For AA < Aniu/e the mixing angle is between 
7r/4 and zero, /5 ~ e^^. Varying AA one goes smoothly from one regime to another. 
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Figure 4: The ratio of tlic integrated rate to the equilibrium concentration of the singlet fermions for 
F = Fq as a function of temperature (in GeV) . The system enters in thermal equilibrium when this 
ratio is equal to one. Left panel: M = 0.14 GeV, right panel: M = 4 GeV. 



We define the temperature T+ at which the singlet fermion enters in thermal equilibrium 
from the equation 

which tells that at T = r+ the number of created particles is equal to the equilibrium one Yeq. 
If <S'+(T) > 1, the initial deviations from thermal equilibrium are damped as exp(— 5+(T)). 
The behaviour of the integrated rate S^(T) as a function of temperature is shown in Fig. [H 
In full analogy, the temperature T_ at which the singlet fermions go out of thermal equi- 
librium is determined by 

S-(T-)^^r(T'-^)'^ = l. (5.35) 
^ ' Y„{T-)h \ dTj T ^ ' 

We show the temperatures T+, T_ and the temperature at which the rate is maximal in Figs. 
[51 The temperature T+ (given roughly by T+ ~ TEw{0.02KM/e)^ at 100 GeV < T < 300 
GeV) is below the sphaleron freeze-out temperature Tew ~ 175 GeV (we take Mh = 200 
GeV) for 

e > 0.02k^ . (5.36) 

Looking at Figs. [3l S] and [5] one can see that thermal equilibrium exists for the range of 
temperatures T_ < T < T-^. As a numerical example let us take the minimal possible mass 
M = and minimal value of Yukawa coupling, F'^ ~ 10~^^ (see Appendix B). It corresponds 
to the choice e = 1 and leads to = ^33- Then the solutions to eqns. (j5.34l5.35p are T+ ~ 15 
GeV and T_ ~ 2 GeV telling that the system is in thermal equilibrium for temperatures 2 
GeV < r < 15 GeV. Asymptotically, the integrated rate approaches 5+(T_) ~ 58. 

Since the Yukawa coupling chosen for this example is the minimal possible one, we reach 
an important conclusion that the reactions associated with the Yukawa coupling ha2 were 
certainly in thermal equilibrium during some stage of the universe expansion. Moreover, 
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Figure 5: The temperatures (in GeV) T+ (upper curves), T_ (lower curves) and the peak rate 
temperature (central curves) as a function of singlet fermion mass (in GeV). Upper panels: normal 
hierarchy, lower panels: inverted hierarchy. Left panels: e = 1, right panels: e = 0.1. 



in the second and third scenarios for the mass difference between singlet fermions with 
AA > Amjy/e the same conclusion is valid for all elements of the density matrix 5pj^ due to 
rapid oscillations between N2 and states. 

The case when e <C 1 is somewhat more delicate. At first sight one may choose e in such a 
way that the rate Tas found in (I5.3ip is always smaller than the rate of the universe expansion. 
And, indeed, the part of it, proportional to Yukawa coupling haz (see (I5.3ip ) is smaller than 
the Hubble rate H for all temperatures if 

6 < 3.4x10-3^ I ^ Normal hierarchy ^^ ^^^ 

M I 0.36 Inverted hierarchy 

Since the mass of singlet fermion is bounded from below by the pion mass, we get that this 
can only happen at e < 2.4 x 10^^. At the same time, the mixing angle (3 gets large at 
T Tp. So, if Tp > T^, the system equilibrates even if (|5.37[) is satisfied. This does not 
happen only if the zero temperature mass difference of singlet fermions is very small, as in 
Scenario Ila, 

< 8 X 10"^- . (5.38) 

To summarise, for any values of parameters, consistent with observed pattern of neutrino 
oscillations, with the exception of the Scenario Ila and for M <^ Mw, the CP-even de- 
viations from thermal equilibrium are damped in some temperature interval [T+,r_] below 
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the electroweak scale. The ratio of the peak rate for equihbration of any element of 6pj^ to 
the Hubble rate is at least 58. If the Scenario Ila is realized, and relations (|5.37l5.38p are 
satisfied, deviations from thermal equilibrimn in CP-even perturbations are substantial for 
all temperatures. Moreover, in any scenario for singlet fermion mass difference, the coher- 
ence in N2 ^ N3 oscillations is lost in the temperature interval [T+,T_]. Thus, the lepton 
asymmetry generation may occur either above r+ or below T_ (see Sec. [7] for details). 

As we discussed, the CP-even deviations are important for generation of the lepton asym- 
metry. The produced asymmetry must not be diluted by reactions that can change it. Thus, 
we consider the CP-odd deviations from thermal equilibrium in the next subsection in order 
to understand whether the asymmetry that was generated before T ~ T-^ or below T_ can 
survive the subsequent evolution. 



6. CP-odd deviations from thermal equilibrium 

The CP-odd deviations from thermal equilibrium are described by eq. ()4.8p . Having found 
the matrices Hint and Pjv in the previous subsection we still should compute six 2x2 matrices 
r^, and 3 rates Pg. They are coming from imaginary parts of the diagrams shown in Fig. 
[6] and have the following structure (we integrated the rates over momenta but are keeping 
the same notations): 



fa ^ J_( \ha2\^R{T,M)-\h^3\^RMiT,M) hl2K^R{T,M) \ 

~ hl^K2R{T,M) \K^\^R{T,M) - \K2? Rm{T, M) ) ' 

n (6.1) 

Tl ~ l^{\h^^\-^ + \h^^\'^){R{T,M) + RM{T,M)) . 
^0 

The minus signs in eqns. ()6.1I) in front of mass corrections come about since the corresponding 
terms in (|4.3l4.4l4.5p are proportional to the chemical potentials pa (notice the change of 
direction of the fermionic line in Fig. [6|) . 

6.1. Approximate conservation laws and damping rates 

The structure of (j6.ip is almost uniquely fixed by the field-theoretical consideration presented 
below. Indeed, the CP-odd deviations from thermal equilibrium can be considered as aver- 
age values of the densities of fermionic currents, which may be exactly conserved for some 
particular choice of the parameters of the z^MSM. 

In the limit when all Yukawa couplings and Major ana masses of singlet fermions are equal 
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Figure 6: The propagator-type diagrams for computation of the damping rates. The Higgs hne can 
be cut and replaced by 'y(T)^; the active neutrino propagator contains one-foop corrections. The 
incoming (outcoming) fermions correspond to two arrows entering (exiting) the vertex. Outcoming 
antifcrmion corresponds to arrows in opposite directions, = iV2 + A^|. 



to zero the z^MSM has five conserved leptonic numbers: 

= Id^xjo , (6.2) 
where a = 1, ...,5. Three of the currents are related to the active leptonic flavours, 

= [L^^L^ + ^„7^^a] , (6.3) 

where Ea are the right charged leptons. The other two conserved currents count the asym- 
metries in singlet fermions and N^, 

= N2l^N2 , A^3^ = Ns-f^Ns . (6.4) 

When the Yukawa couplings and Majorana masses are switched on, none of these numbers 
are conserved any more. 

To make the discussion more transparent, consider the following combinations of the cur- 
rents introduced above: ^ 

a=l 

and ^ 

4 = J^ = Y.J' + ^2+N!^ ■ (6.6) 

a=l 
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The first current Jj^ (total leptonic number) corresponds precisely to the leptonic number 
symmetry defined in [7] which is exact in the limit has 0, AMjj — > 0, whereas the second 
current (it can be called total fermionic number) is conserved when all Majorana neutrino 
masses are put to zero. What concerns the currents for a given a, they are conserved in 
the limit haS — > 0, ha2 0. 

Now, if some combination of the currents introduced above is exactly conserved, the equa- 
tions (14. 8p with zero source terms must have a time-independent solution for any choice of 
initial conditions. As an example consider first the limit M ^ 0, e ^ 0. In this case the 
current is exactly conserved, and we must have 

d_ 

Iti 

for any (5/?_ and /^q. This leads to 



Tr(5p_ 







(6.7) 



TN = ^fl and Trf^ = r2 for M = 0. (6.J 
a 

In the another limit M 7^ 0, e — > it is the current J4 which is exactly conserved and 



d_ 
di 



Trr3(5p_ +^1^-0 



{6: 



for any Sp- and pa (here T3 is the Pauli matrix). This gives 



-p22 

J- N 



a22 
L ) 



-.33 

■ N 



a33 
L 1 
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= and TrrgT^ = for e = 0. (6.10) 



In more general terms, the consistency condition can be formulated as follows. Rewrite 



eqn. ()4.8p with = = in the form 



dz 
'dt 



Dz , 



(6.11) 



where z is a vector with 7 components, z = {5p_ , 5p_ ,Sp_ , 5p_ , p-a) and D is the 7 x 7 matrix 
constructed from Tg, T^, Tg and H. The time-independent solution appears when D has 
a zero eigenvalue. Then, we must have det D = for the following choices of parameters, 
corresponding to the conservation of the 5 currents introduced above: hu = for / = 2, 3, 
corresponding to conservation of the leptonic number of the first generation (and similar 
relations for the second and third generation), ha3 = 0, corresponding to conservation of the 
current J4 (and an equivalent relation for ^ -^3)5 and M = 0, leading to conservation of 
J5 . One can check that eq. (j6.1|) indeed satisfies these requirements. 

It is instructive to find the damping rates in the limit M — > 0, e ^ 0. In this case the 
matrix D has two zero eigenvalues corresponding to the conservation of currents J4 and J5 , 
2 complex eigenvalues 

llR(T,M)±i{E2-Es) (6.12) 
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corresponding, as in the case of CP-even perturbations, to the off-diagonal elements of the 
density matrix 5/9_, and three eigenvalues related to the damping rates of three different 
leptonic flavours, 

7^ = ^i?(T,M), (6.13) 
where the roots of the cubic equation 



^2 3 / _E>^\ IMAM^ 



x^ + 2x^ + ^i^l-^^jx+ -y/- =0. (6.14) 

If, for example, /ie2 ^ ^fi2> ^e2 ^ ^t2 then the smallest root of eq. (|6.14|) is approximately 
given by 5/ig2/4-F^. From (j6.12|) we can see that the coherence in CP-odd perturbations is 
lost at the same time as it is in CP-even perturbations. As for the damping rates of active 
flavours, with the use of constraints ()B.2IB.3P (see Appendix B) one finds that the integrated 
rates corresponding to 7^ 

1 . , dT 



are at least 



Si ~ 8.2/e, S2 ^ 50/e, S3 ~ 156/e , Normal hierarchy , (6.16) 

Si ~ 32/e, 52 ~ 22/e, S3 ~ 122/e , Inverted hierarchy , (6.17) 

where the smallest number in (j6.16p corresponds to the asymmetry in the electronic flavour. 
Eq. (|6.17p shows that if the hierarchy is inverted, all the rates exceed the rate of the universe 
expansion by a factor of at least 22 (corresponding to the damping of asymmetry which existed 
before the equilibrium period by a factor smaller than ~ e" -22 rs. 3 X 10"^°). For the case of 
the normal hierarchy eq. (j6.16p shows that the damping is at least ~ e"^'^ ~ 3 x 10~^. This 
leads to the conclusion that the reactions which change leptonic numbers in each generation 
were certainly in thermal equilibrium during some time below the electroweak scale which is 
good enough to dilute the lepton asymmetry below the level required for resonant production 
of dark matter. At this point the z/MSM is very different from the Standard Model, where 
leptonic numbers are conserved (up to electroweak anomaly). 

6.2. Protection of lepton asymmetries 

The fact that the flavour changing reactions were in thermal equilibrium during some period 
of the universe expansion below the electroweak scale would at first sight mean that no 
(large) asymmetry in active leptonic flavours can exist at small temperatures. However, this 
conclusion is not necessarily true since some combination of asymmetries in active and sterile 
flavours may be protected from erasure due to the existence of approximate conservation laws 
of currents J4 and J5 . The only certain thing for the moment is that the low temperature 
remnants of high-temperature leptonic asymmetries in active neutrinos are flavour-blind, i.e. 
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/^e ~ /^^t — /Ur = /fi- This fact allows to simplify the further analysis replacing the system of 
equations ()4.8p with zero sources by 

.d5p- ITT s: 1 ^ 



= [H,5p^]--{TN,5p^] + ti^Y.^% ^ 



dp 1 



3 . 



(6.18) 



To consider the possibility of protection of lepton asymmetry we start from the Scenario 
I for the mass difference of singlet fermions. Then for bp ^ the Hamiltonian Hint can 
be diagonalized simultaneously with the damping rates in eq. (I6.18p . and one finds that for 
small e and M the rates 74 and 75 are 

74 ~ r^R{T,M) , (6.19) 



75 - --^Rm{T,M), (6.20) 



4^2 

where is defined in (|5.32|) . In comparison with 71,2,3; the rate 74 is suppressed by whereas 
the rate 75 is suppressed by M^. For bp < the mass matrix Hint is not proportional to 
Fat, J2a ^'^d J2a longer, leading to the rate 

74^74 + sin2/3r22 , (6.21) 

where r22 is defined in eq. (j5.3ip and the angle /3 in eqns. (j5. 915. lip . At the same time, the 
rate 75 is not changed. Now, repeating the considerations of the previous section one finds 
that is protected from erasure only if inequalities (|5.37l5.38p are satisfied simultaneously, 
i.e. only for Scenario Ila. 

Another leptonic charge which can be protected from erasure by the processes with lepton 
number non-conservation is J^. If max (^) < 1 then the density matrix at low temperatures 
has the form 

Peg = exp -^sQsj = exp -;U5(L + (52 + Q3)j , (6.22) 

where is the chemical potential corresponding to the effectively conserved charge = 
fd^xJ^. In this case the previously generated asymmetry in Q5 survives, and the fact that 
Q5 contains the currents corresponding to active flavours ensures non-zero asymmetry in 
lepton number, which is essential for the resonant production of dark matter sterile neutrinos. 
Independently of the choice of parameters, the chemical potentials for are the same 
as those of the active fermions, which is the consequence of the fact that the transitions 
La — > A^2,3 are in thermal equilibrium. Exactly the same conclusions are valid for the second 
and third scenarios for the fermionic mass difference. 

The region of the parameters in which the asymmetry in Q5 is protected can be found 
from the condition that the peak value of j^/{2H) does not exceed 1. We plot this region in 
Fig. El 
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Figure 7: The region of parameters in [e (vertical axis), Af/GeV] plane for which the low temperature 
lepton asymmetry is "protected" from erasure for normal (left panel) and inverted (right panel) 
hierarchies of neutrino masses. The upper curve corresponds to the damping factor e~^, the lower 
curve to 0.002, and the middle one to 0.1. 

To summarize, the existence of a lepton asymmetry at small temperatures ~ 100 MeV is 
only possible in the following situations: 

(i) The asymmetry is produced below the temperature T_, when the processes that damp 
the CP-even and CP-odd deviations go off thermal equilibrium. 

(ii) The asymmetry in is produced above r+ and the i^MSM parameters lie in the range 
shown in Fig. [7] ensuring that it is not erased later on. 

(iii) The asymmetry in Q4 is produced above T4. and the z^MSM parameters lie in the range 
(I5.37I5.38|) . 

In the next section we will add to the analysis an input from the dynamics of lepton 
asymmetry generation which will allow to choose between these possibilities and to add 
further constraints. 

7. Lepton asymmetry generation and constraints on masses and couplings 
of singlet fermions 

To find the leptonic asymmetry one should solve equations (j4.3|4.4l4.5p with zero initial 
conditions for chemical potentials and for the elements of the density matrices of singlet 
fermions. Due to the fact that the number of equations and different time scales is large (the 
equation count for real variables is as follows: 4 for p, 4 for p, and 3 for /Iq,) this cannot 
be done analytically. Nevertheless, the behaviour of the system can be understood on the 
qualitative level with the results of Sec. [5] and Sec. [6] and a number of quantitative estimates 
can be made. 

Let us start from the small time behaviour of the system, when all reactions involving 
singlet fermions are out of thermal equilibrium, so that the largest exponential in (jS.SOp . 
T22t is smaller than 1. This regime was considered in [2] for Scenario III, assuming that 
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the number of oscillations of singlet fermions, defined below in eq. (j7.14p . at the time of 
electroweak cross-over, is much larger than one. We will generalize this analysis to a more 
general case, accounting for the electroweak symmetry breaking effects and considering also 
the time so short that 5M{T)t < 1. 

For these purposes it is convenient to transform the system in a form that does not contain 
the term responsible for oscillations between the two singlet fermion flavours, [H, pi\f], see [2]. 
This can be done by introducing p^- related to p in the following way: 

PN = U{t)E{t)pNE\t)U\t), E{t) = e^p(^-i dt'AE{t')^ , (7.1) 

where the matrix U{t) converts the Hamiltonian Hint to the diagonal matrix AE{t) defined 
ineq. §28^, 

H,nt = U{t)AE{t)U\t) . (7.2) 

Then the equation for pat is 

= pn] - \{TNurPN - /o^"} + ipatU , (7.3) 

where 



H = —E^(u^U-U^u')e 



2i 

Tmu = E^U^TnUE , (7.4) 
f^^ = E^U^t%UE . 

Exactly the same procedure applies for the equation describing the antiparticles. 

As was explained in [2], the set of equations (j4.3l4.4l4.5p can be solved perturbatively for 
the case when all damping rates (symbolically F) are small enough. Ft ^ 1. This is done in 
the following way: rewrite the differential equations (j7.3p in the integral way, e.g. 



PN 



\ I dt'(right hand side of eq.dZS])) (7.5) 



and then solve them iteratively. Then asymmetries in leptonic numbers (chemical potentials 
//q) are given by 









P-a = 


f dt' 


/ dt"Tr 


Jo 


Jo 
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f dt' 


/ dt"Ti: 




Jo 


lo 



where 



{tl{t')V{t',t")TN{t")) V^{t',t") 
{Tf{t')V{t',t")T*^{t"))v\t',t")\ , (7.6) 



V{t',t") = U{t')E{t')E\t")U\t") , 

V{t',t") = U{t'YE{t')E\t")U^{t") (7.7) 
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(T corresponds to the transposed matrix). Equation ()7.6p can be simplified, 



^^{t) = 4 /" dt' [ dt' 



.11 



Im (t/t(tOf£(tO^(tO)i2(C/^(t")r^(i")f^(i"))2i X 



JO Jo 




(7.8) 



As usual, the asymmetry contains a product of two imaginary parts. The first multiplier in 
(j7.8|) is associated with the CP-breaking complex phases in the Yukawa couplings, whereas 
the second corresponds to the oscillations between two singlet flavours. As was shown in 
[2], in the second order of perturbation theory and neglecting mass corrections O (y^)) the 
total leptonic asymmetry is zero, /^a = 0- It appears in the third order only, leading to an 
extra suppression of the order of Ft. In this work we will not go beyond the second order 
of perturbation theory, and account for extra suppression by multiplying the results by Tt. 
Note that the resonant production of dark matter sterile neutrinos occurs even if total lepton 
asymmetry is zero but individual flavour asymmetries are large enough |161 [2] . 

What happens for large times? For definiteness, suppose that the smallest damping rate for 
CP-even deviations from thermal equilibrium is F33. Then, for T^.^t ^ 1, the CP-even fluc- 
tuations thermalize, 5/9+ ^ 1, and the source terms in ()4.8|) completely disappear, meaning 
that the production of lepton asymmetry stops. 

In fact, an even stronger statement is true, namely that there is no generation of lepton 
asymmetry for T22t ^ 1- Indeed, in this regime all but one element of 5pj^ are exponentially 
damped: the oscillatory off-diagonal part of the CP-even density matrix disappears at (F22 + 
F33)t/2 > 1, and one of the diagonal elements at T22t > 1. So, in the mass basis 



For this type of deviation from thermal equilibrium the source terms in the equation for 
chemical potentials ;Uq,, eq. ()4.10p vanish, Sa = 0. The same is true for diagonal elements of 
the density matrix accounting for asymmetries in singlet fermions. In other words, the 
leptogenesis ceases to work when coherence in oscillations of singlet fermions is lost, which 
happens when one of them is thermalized. The same conclusion is reached if the damping of 
coherent oscillations is inserted "by hands" into eq. (j7.8p . 

To conclude, we expect that the asymmetry is maximal at tcoh ~ 2/(Fi +F2). For t > tf-oh 
the production of the asymmetry is switched off, and the asymmetries in different quantum 
numbers decay with the rates found in Section [HI 

Let us estimate the maximal possible asymmetry which can be created at t ~ tco/u corre- 
sponding to the temperature at which A''2 equilibrates, T ~ r+. 

In the Scenario III for the singlet fermion mass differences the matrix U{t) depends on 
time slowly. Indeed, when the tree level mass difference is much larger than the Higgs induced 




(7.9) 
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mass, the matrix U{t) corresponds to the rotation by 7r/4, and U ~ Ami,/ AX. Therefore, 
the asymmetries at time t are of the order of 

^a{t) ^ 5cP-^m, (7.10) 

^0 



where 

expU / dt"'{E2{t"')-E^{t"')) 



$(t) = dt' f dt"R{T',M)R{T",M) x Im 
JO Jo 




(7.11) 



where the temperatures T', T" correspond to the times t' , t" , and 6cp is defined in (|3.19p . 
For e ~ 1, (5cp can be of the order of 1. 

In the Scenario I for T > Tp the mass difference is determined by the vev of the Higgs 
field only. Therefore, the temperature dependence of the matrix U (t) can also be factored 
out up to mass corrections M'^/T'^, so that U ~ M'^/T'^. However, the asymmetry in fi^ is 
suppressed in comparison with eq. (j7.10p by a factor (at small e) 



|/lt/l|22 / \Kmatm 

since in the limit e — > the matrices Hint and F^v can be simultaneously diagonalized (cf. 
eqns. ()5.22p and (j5.2p ). so that off-diagonal elements appearing in (j7.8p are suppressed either 
by a factor Sj or by a mass to temperature ratio M"^ /T"^. A similar factor appears in the 
Scenario II for e <C 1. 

For the generic case of Scenarios II the phase factor cannot be factored out and the 
equations are more complicated. We expect, however, that the discussion below has a general 
character, at least on the qualitative level. 

It is instructive to find the behaviour of ^{t) in different limits. For this end we will 
assume that the rate -R(T, M) can be approximately represented as R{T, M) = AT~'^, where 
n is some number. For example, for temperatures above the peak of production of singlet 
fermions n ~ 4, at T > 100 GeV n ~ 1, whereas at temperatures below the peak n ~ —5. 
The exponential in (j7.10p can be written as 

dt"'{E2{t"') - E^{t"')) = x{T') - x{T") , (7.13) 



Jt" 



where 

= [ at [^^) . o.i,"'"^^^"" . (7,14) 

and (...) is the thermal average. The physical meaning of the parameter x{T) is that x{T)/2tt 
gives the number of oscillations between singlet fermions from the end of inflation till the 
temperature T. Then one easily finds: 

$(i) = const / dzizf^ / dz2sm.{zi — Z2)z2 
Jo Jo 
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Figure 8: The behaviour of functions n = 4 (left) and F^{x), n = —5 (right) counting the 

number of singlet fermion oscillations near the temperatures T+ and r_ . 



where -^+(2;) in limiting cases is given by 
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F+{x) 



(n+2)(n+5)(2n+7)' 
3 1 
2?i+l X 



X <C1 

X > 1 



(7.16) 



which is valid for n > — 1, true for any temperatures T > Tmax, where T^ax is the temperature 
at which the rate of production is maximal. The plot of the function i*+(x) for n = 4 is 
shown in Fig. [8l 

Due to the very steep dependence of <I>(t) on the temperature the baryon asymmetry, 
produced at T ~ T^w can be much smaller than the lepton asymmetry, created at T ~ r+. 
Indeed, for n = 4 and for x > 1 one gets that fia oc l/T^^, so that a drop of the temperature 
by just a factor of 2, increases the asymmetry by a factor of 3 x 10^. Including an extra factor 
Tt ~ R{T)/3H, accounting for the fact that baryon asymmetry is produced in third order of 
perturbation theory |2J amplify the difference even further. 

Let us estimate the maximal possible asymmetry which can be produced at T+. For this 
end suppose that the number of oscillations maximizes the function (^pmax ^ 0.076 at 
X ~ 3.8) and that CP-violation is maximal. Clearly, A cannot be larger than Amax = 
4/(9 X 2 + 4) = 2/11, where 4 is the total number of spin-states of and 9 is the number 
of spin-states of three leptonic generations. Thus, 



A ~ Ar 



eF+(x(r+)) 

' pmax 



(7.17) 



where the factor e accounts for the fact that CP-violation goes away in the limit e — > 0. 
Similar estimates apply for asymmetries in the other quantum numbers defined in Section 



m 



5Q4 



e2F+(x(r+)) 



FV 



(7.18) 



An extra factor e appears since the rate of creation or destruction of is suppressed by 
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in comparison with the rates changing As for the asymmetry in Q5, one gets 

SQ,:.A^..^^^^'J±^^, (7.19) 

where the second factor takes into account that the processes with the change of are 
suppressed in comparison with L <-> A^2 transitions. 

The asymmetries in different quantum numbers generated at T ~ T+ are reduced later 
with the rates determined in Section [6l 

7.1. Constraints on singlet fermions from baryon asymmetry 

The estimates of the leptonic asymmetry presented above allow to find constraints on the 
masses and couplings of the singlet fermions from the requirement that the produced lepton 
asymmetry is large enough to make baryon asymmetry at the freezing point of sphaleron 
processes. 

Consider first Scenarios I, II for the singlet fermion mass difference. If r+ > Tewi the 
asymmetry generation in this case occurs in the resonant regime as the number of oscillations 
at temperature r+ does not depend on parameters and is of the order of one, 

AT.) . 12^ . 

We present in Fig. [9] the region of the parameter space in which the baryon asymmetry 
()7.17p . damped by a factor ex.p{—S+{TEw)) can exceed the observed value for normal and 
inverted hierarchies. We take the sphaleron freeze-out temperature to be 175 GeV, corre- 
sponding to the Higgs mass 200 GeV [62j and account for a suppression factor Sj defined in 

eq. irmn . 

The asymmetry related to the charge (54, eq. (j7.18p can exceed the observed baryon 
asymmetry for e > 10~^ for normal hierarchy and for e > 10^^ for a wide range of the 
singlet fermion masses, including M > Mw. The fact that the baryon asymmetry generation 
is also possible for masses so large was missed in [2] and is due to the fact that the charge 
Q4 is protected from erasure for small e, whatever the value of M is. 

In Fig. [in] we show the region of the parameter space where the asymmetry in can 
exceed the observed value for the case of the normal hierarchy. The parameter e is bounded 
from below by e ~ 7 x 10~^, and the mass from above by M ~ 100 GeV. These results refine 
the estimates presented in [3. 

In the Scenario III the leptogenesis goes off the resonance and the available parameter 
space decreases. Comparing eq. (|7.17p with observed baryon asymmetry one can put an upper 
bound on the mass difference difference of singlet fermions, 5M/M < 4 x 10^^K'^(Af/GeV), 
valid if < Tewi M < 50 GeV and e 1. If lies in the symmetric phase of the 
electroweak theory, > 250 GeV, a constraint from [2j, 5M/M < 6 x 10~®(M/GeV)t 
should be used. 
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Figure 9: The region of the parameter space in [e (vertical axis), M/GeV] plane in which the asym- 
metry defined in eq. (|7.17p . and then reduced due to damping, can be consistent with observations. 
The lower line corresponds to asymmetry A = 6.6 x 10~^ (corresponding to observed baryonic asym- 
metry), the middle one to A = 6.6 x 10~^ and the upper line in left panel to A = 6.6 x 10^'^. Left 
panel - normal hierarchy; right panel - inverted hierarchy. 

7.2. Low temperature lepton asymmetry 

Let us find now the region of parameters which can lead potentially to the generation of a 
large lepton asymmetry (AL/L > 2x 10~^, as required by observational constraints, discussed 
in [39]). Clearly, the constraints coming from baryon asymmetry are much weaker than those 
related to the large lepton asymmetry at lower temperatures. As we have already discussed, 
the asymmetry can be generated somewhat above r+ or below r_. 

We start from T ~ T+. Out of five different leptonic numbers discussed in Section [6] only 
two can survive the subsequent evolution. These are the asymmetry in in the Scenario 
Ila, provided e is small enough and in Q5 which is protected if the mass of singlet fermion is 
small enough, see fig. [71 

As we saw the number of oscillations at T+ plays an essential role in the determination of 
the asymmetry. So, we present in Fig. [TTjthe quantity x(T+) for the Scenario I of singlet 
fermion mass difference for e = 1 for the case of normal and inverted hierarchies. For the 
generic choice of parameters for the Scenario II the number of oscillations is of the same 
order. However, by tuning the Majorana mass difference to the Higgs induced mass difference 
it can be made much smaller (for e = 1, see eq. ()5.4p ) than the numbers appearing in Fig. [TTl 
For the Scenario III the number of oscillations is much larger than that in the Scenario I 
(by a factor 5M/matm if the comparison is with normal hierarchy case). 

Consider now the lepton asymmetry in Scenarios LIIL 

Scenario I. The only possibility is to have an asymmetry in Q5. Inserting different 
rates in (I7.19|) we get for e ~ 1 the asymmetries plotted in Fig. [12] (assuming that the 
number of oscillations maximizes the asymmetry) . For the normal hierarchy the asymmetry 
does not exceed 2 x 10~^ and thus is smaller than the minimal required number (2 x 10~^) 
at least by a factor of 10. For the inverted hierarchy the maximal asymmetry is about 
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Figure 10: The region of the parameter space in [e (vertical axis), A//GeV] plane in which the 
asymmetry defined in eq. (|7.19p and reduced later due to damping discussed in Section [SI can be 
consistent with observations. The upper line corresponds to asymmetry A — 6.6 x 10"^, the middle 
one to A = 6.6 x 10~^ and the lower line to A = 6.6 x 10~^. We took the normal hierarchy case. 

1 X 10~^, a factor of 20 smaller than required. Though there are no orders of magnitude 
differences between potentially produced asymmetries and the required one, the conclusion 
that Scenario I cannot lead to necessary lepton asymmetry is robust. Indeed, in all estimates 
the CP-violating affects were assumed to be maximal, and other uncertainties were pushed in 
the direction which can only increase the asymmetry (for example, accounting for the number 
of oscillations will reduce the asymmetry for the case of normal hierarchy by a factor of 20). 

In the Scenario II for generic choice of parameters the results for stay the same as 
in the previous case. In other words, no sufficient asymmetry in can be produced at 
for this case. Potentially, in the Scenario Ila the leptonic charge Q4 can survive. However, 
this can only happen if e < 2.4 x 10^^. For e so small the maximal asymmetry in Q4 cannot 
exceed Amax^^ — 10~^, too small to have any effect on dark matter production. Now, if the 
Scenario III for singlet fermion mass difference is realized, the asymmetry gets reduced by 
a factor niatm/^M <C 1 in comparison with Scenario I. Since no large asymmetry can be 
produced in the Scenario I, Scenario III can be discarded as well. 

To summarise, no generation of large lepton asymmetry at T ~ T-f, which can survive till 
small temperatures, is possible. 

Consider now a possibility of large lepton asymmetry generation at lower temperatures, 
T ~ r_. The oscillations of singlet fermions re-enter into coherence regime at T ~ T_, 
corresponding to t-. Then, one can simply change the region of integration in (|7.8p : 



accounting for the fact that at t < t_ the oscillations were exponentially damped. Cor- 




(7.20) 
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Figure 11: The number of oscillations (vertical axis) of singlet fermions at temperature T+ for the 
Scenario I as a function of the fermion mass (in GeV) for normal (upper red curve) and inverted (lower 
blue curve) hierarchies. We took e = 1. 
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Figure 12: Maximal possible lepton asymmetry (vertical axis) generated at T+ and survived till T_ 
for the Scenario I for normal (left) and inverted (right) hierarchies as a function of singlet lepton mass 
(in GeV). We took e = 1. 



respondingly, the limits of integration in the phase factor <I>(T) defined by (jT.lip must be 
changed. We get: 



$(i) = const X Im / dzie''' zi""^^^^ / dzae-^'^zf"^^^^ 

Jx(T) Jx(T) 



lx{T) ^ Jx{T) 



2 



= (^)'^-(^(^))' (7-21) 

where the plot of the function F^[x) for n = —5 is shown in Fig. [8l It reaches the maximal 
value F™"^ = 0.167 at x = 0.47. 

In limiting cases the function F-{x) is given by 

-(„+2)(„+5)(2n+7)^ ' ""^^ (7.22) 
~2^x ' ^ > 1 
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which is vahd for n < —5. The case of n = —5 requires a special treatment, leading to the 
asymptotic value F„(x) = — xlog(e7x) for x ^ 1, where e = 2.718... and 7 = 0.577... is the 
Euler constant. 

To estimate the leptonic asymmetry generated at this time one can write 

eF_(x(r_))5n 

A ~ Amax pmax ~ ' (7.23) 

where the factor 6n /riu accounts for deviation of the sterile neutrino concentration from the 
equilibrium one at temperatures close to but below T_. If A^2,3 decouple from the plasma being 
relativistic, M < T_, the deviation of their concentration from equilibrium is suppressed by 
the factor 

5n/nu ^ \negiT,M)/neqiT,0) - 1| ~ 0.2M'^/tI . (7.24) 
If the decoupling occurs when T- < M, the corresponding factor is 



Since all reactions which change different leptonic numbers are out of equilibrium at temper- 
atures below T_, the asymmetries ()7.23p stay intact. 

Let us estimate the number of oscillations at t ~ t_. Suppose first that T_ ^ M so 
that the high temperature approximation can be used. Then with the use of eq. (j5.27p the 
temperature T_ is given by 

T ( !M ^ (7 26) 

\KBGlmatmMo) ^ ' ' 

leading to 

x{t-) ^ ^L^{GFMofmatm6M . (7.27) 
e 

The asymmetry is maximal if the number of oscillations is minimal. So, to get the maximal 
asymmetry we should take e = 1 and the minimal 5M. For the Scenario I this corresponds 
to the inverted hierarchy of neutrino masses and to x{T-) ~ 3.6 x 10'^. So, the asymmetry 
cannot exceed 10~^, with the actual number being smaller as one has to account for the factor 
Sn/riy < 1 and for extra suppression from CP-breaking phases. For the normal hierarchy of 
neutrino masses the number of oscillations is larger by a factor of ~ 50, and for the Scenario 
III for the singlet fermion mass difference it is even higher. We conclude, therefore, that large 
lepton asymmetry, interesting for dark matter production, cannot be generated at T ~ r_ 
for Scenarios I and III, at least if M ^ T_. 

Let us find the critical singlet fermion mass where the relativistic approximation used above 
is not valid. Since the typical momentum of a fermion in the plasma is {p) ~ 3T, we require 
3T ~ M and find that the singlet fermions decouple being non-relativistic if 

/ 27c V'^ 



38 



giving Merit — 30 GeV for e = 1. We will demonstrate now that the lepton asymmetry is also 
very small if the singlet leptons decouple in the non-relativistic regime (again Scenarios I 
and III are considered). 

At large singlet fermion masses one can neglect the influence of the medium and consider the 
processes involving A''2,3 as if they were in the vacuum. The fastest reactions at temperatures 
T < M\y are the decays Z — > uN and W IN (with the rate Fy), and decays or inverse 
decays of N to all possible leptonic or semi-leptonic channels (rate Tjsr). The rates of inverse 
W and Z decays, responsible for thermalisation, can be approximated as (at M < M^r) 




M2 




3/2 



(7.29) 

where T]y_^ii^ ~ 0.7 GeV and Tz^pu — 0.5 GeV are the widths of the intermediate vector 
bosons, ny = 3, and 9q ~ Kmatm/C^^M). 

The rate of inverse decays of N is of the order 

where A is proportional to the number of open channels for decays, ^ ~ 10 if M > 10 GeV 
|llj . The temperature at which the oscillations of N start to be coherent can be determined 
from the condition H = Ti\f + Ty , and the lepton asymmetry from the relations (|7.23|7.25p . 
The results for the temperature T_, the number of oscillations and the lepton asymmetry are 
shown in Fig. [T3l Note that for the non-relativistic case the number of oscillations is given 

(7.31) 

One can see that the asymmetry never exceeds 2 x 10~^ (inverted hierarchy) and 4 x 10^^ 
(normal hierarchy), which is well below the threshold for the resonant production of dark 
matter. We conclude, therefore, that no substantial asymmetry generation can occur after 
singlet fermions decouple in Scenarios I and III. The same conclusion is valid for the 
Scenario II for a generic choice of parameters. 

On the other hand, for a special case of Scenario Ila, when the sterile fermion mass 
difference is much smaller than the active neutrino mass difference, the asymmetry production 
enters into resonance and the generation of large lepton asymmetries AL/L > 2 x 10~^ at 
T_ becomes possible for a variety of masses and couplings of singlet fermions. With the use 
of eqns. ()5.7|7.14|7.26|) one finds that if 5M(0) = 0, the number of oscillations at T = T_ 
does not depend on M and e (for e ^ 1) and is given by 

p\'^ 8 7iT 
TTaw 



x{T^)^0.15(n ^_^(2 + cos2M^10. (7.32) 
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Figure 13: The temperature T_ in GeV (upper left), the number of oscillations (larger for normal 
hierarchy) of singlet fermions at temperature T_ (upper right) and the maximal lepton asymmetry 
(smaller for normal hierarchy) generated at T_ as a function of the fermion mass (in GeV) for normal 
(red curve) and inverted (blue curve) hierarchies. We took e = 1 and the Scenario I. 



In other words, we are close to the resonance and a large asymmetry can be produced. In 
Fig. [13] we present the part of the parameter-space where the asymmetry may exceed the 
critical value. 

In fact, yet another mechanism for late leptogenesis is possible in the Scenario Ila with 
the "tuned" mass difference. If 3r_ > M, the singlet fermions decouple from the plasma 
being relativistic. Later, they decay with lepton number non-conservation and CP-violation 
and, if they are degenerate enough, they will produce large lepton asymmetriej§. 

Let us estimate the value of the lepton asymmetry which can be created in decays of 
Since the Yukawa couplings are very small, the main contribution to CP asymmetry comes 
from the mixing between A''2 and N^, as shown in Fig. [T5j An estimate for asymmetry reads 



A ~ A. 



max 



\ 6M ; 



(7.33) 



_rivMo_ 

where the the first term (e) comes from CP- violation, the second term describes the resonance 



*The fact that CP-violation is greatly enhanced in the decays of degenerate particles is well known from 
physics. It was first suggested for baryogenesis in [63], discussed in [64] and studied in detail for TeV scale 
Majorana fermions in [651 1661 167] . 
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Figure 14: Left panel: the parameter-space (I and II) which can lead to the lepton asymmetry, 
produced at T = T_ and exceeding 2 x 10"'^. Right panel: the parameter-space (I and II), which 
can lead to the lepton asymmetry, produced in decays of N2.3 and exceeding 2 x 10^^. In the region 
III A^2,3 decay below the temperature 100 MeV and thus do not contribute to resonant production of 
dark matter. In the region I (II) A^2,3 decouple being relativistic (non-relativistic). 
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Figure 15: Diagrams for A'2,3 decay which can lead to large lepton asymmetry below the electroweak 
scale. 



and is valid for 6M > F^r (it should be replaced by 1 in the opposite limit), the third term 
accounts for equilibration of the asymmetry due to inverse decays. It should be replaced 
by 1 if M2 > Fat Mo, i.e. for 



For the case 5M ~ F^r the asymmetry can be large and lead to the resonant production 
of dark matter sterile neutrino, provided A''2,3 decay above the temperature ~ 100 MeV, at 
which A^i are created most effectively. The latter requirement leads to the constraint 



GeV > 'H^inFsJ Uj ■ P''^) 

In Fig. [14] we present the part of the parameter-space where the asymmetry created in A^2,3 
decays may exceed the critical value. 
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To constrain further the parameter-space of the model one should take into account the 
requirement that not only the low temperature lepton asymmetry must be large enough, 
but also that the baryon asymmetry is small. Given the number of CP-phases and other 
parameters we expect that anywhere in the regions shown in Fig. [14] the required hierarchy 
can be achieved by some choice of Yukawa couplings. However, for a generic case, in which 
no cancellation between different CP-violating phases takes place, the region of small singlet 
fermion masses and large e is singled out. 

Indeed, in the Scenario Ila the baryon asymmetry generation occurs in the resonant 
regime at Tj^ > Tew-, leading generally to large baryon asymmetries. The observed small 
baryon asymmetry can be derived moving out of the resonance, i.e. for Tj^ < Tew- The 
number of oscillations at the electroweak temperature Tew — 175 GeV is (for e ^ 1): 



and smaller than one if e is large and M is small. In this regime the baryon asymmetry is 
suppressed by a factor 



which is about 5 x 10~^ for M ~ 2 GeV, e ~ 1 ,k = 1 (upper left corner in Fig. [T^ . 
producing roughly a correct hierarchy between high temperature baryon asymmetry and low 
temperature lepton asymmetry. 

Finally, let us discuss the possibility that large lepton asymmetries Aq ^ 2 x 10^'^ were 
generated well above the electroweak temperature. Is it possible that they were not trans- 
ferred to baryon asymmetry but survived till low temperatures? 

As we have already found, the only leptonic numbers that can survive till low temperatures 
are related to the currents and J^, defined in (|6.5|6.6p . Moreover, the only flavour structure 
of primordial asymmetry which is consistent with small baryon asymmetry is the one in which 
L + AN2 = 0, where L is a lepton number of active fermions, and is the asymmetry in 
a more strongly interacting singlet fermion. Indeed, if L is large it will lead to large baryon 
asymmetry due to sphalerons. If is large, a part of it will be transferred to L and then 
to baryon asymmetry. The amount of going to L is at least 



where we used the minimal possible rate R{T,M) corresponding to M ~ m.^ and e = 1. In 
other words, the only possibility is to have large asymmetry in A^3, AN3 = Aq and assume 
that e ^ 1, suppressing the transitions N3 ^ L. 

Now, four different possibilities can be realised. If the reactions changing and were 
both in thermal equilibrium, no primordial asymmetry will survive. If, on the contrary, none 
of the reactions changing (^4 and Q5 were in thermal equilibrium, a large asymmetry in 




(7.36) 




(7.37) 



AoS+{Tew) > 2 X lO^^Ao > AB 



(7.38) 
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Figure 16: Left panel: Part of the parameter space corresponding to conservation of Q4 and non- 
conservation of Q5. This can only be realised in Scenario Ila in which Tp < Vertical axis: e, 
horizontal axis: mass in GeV. The admitted regions are below the curves. Upper red line - normal 
hierarchy, lower blue line - inverted hierarchy. Right panel: Part of the parameter space corresponding 
to conservation of Q5 and non-conservation of (54. It is required that Tfj > T_. Vertical axis: e, 
horizontal axis: mass in GeV. The admitted region is to the left of the curve. No parameter space is 
allowed for the inverted hierarchy case. 

N3 will not be transferred to an asymmetry in active leptons, and, therefore, no resonant 
production of dark matter sterile neutrinos is possible. So, to get large lepton asymmetry at 
low temperatures one must require that one of charges out of Q4 and must be conserved 
and the other equilibrate. In Fig. [16] we present the parameter-space in which the primordial 
asymmetry in N3 induces a baryon asymmetry smaller than the observed one but leads to 
large low temperature lepton asymmetry. It requires rather small values of the parameter e. 

8. Fine tunings or new symmetries? 

The requirement that the z^MSM produces both baryon asymmetry and dark matter in 
amounts required by observations puts very stringent constraints on the parameters of the 
model. In this section we will discuss whether these constraints, appearing as different fine- 
tunings in the Lagrangian of the i^MSM, can indicate the existence of some hidden approx- 
imate symmetries. These symmetries, if exist, cannot be explained in the framework of the 
z^MSM itself, as this model is based on a renormalizable field theory which may be valid all 
the way up to the Planck scale [12j. At the same time, their presence can give some hints on 
the properties of more fundamental theory, replacing the z^MSM at high energies. 

We start from the relative strength of Yukawa interactions of singlet fermions N2 and N3 . 
A non-trivial constraint on the z^MSM parameters is coming from the requirement that the 
baryon asymmetry at the electroweak scale must be much smaller than the lepton asymmetry 
at small temperatures. It tells that the parameter e should be close to its maximal value, 
e ~ 1. For e that large the lepton number U(l) symmetry, introduced in [7], is strongly 
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broken in the singlet fermion Yukawa sector, but is respected by the Majorana masses of the 
singlet fermions and by charged lepton Yukawas. Therefore, one may wonder if some other 
global symmetry, respected both by the Yukawa couplings and by Majorana masses, may 
exist for the extreme case e = 1. 

As was discussed in [7J, such a symmetry does not exist if both charged and singlet lepton 
Yukawa couplings are taken into account. If, however, charged lepton Yukawas are disre- 
garded, quite a symmetric singlet lepton interaction can be found in the inverted hierarchy 
case. Indeed, for the case mi = m2, 623 = 7r/4, and ^13 = the fields L2 and L3 defined in 
(|3.9p are the orthogonal mixtures of different leptonic flavours. Thus, for e = 1 the Yukawa 
part of Lagrangian (j3.3p is symmetric with respect to the non-Abelian flavour group SU(2) 
(broken, of course, by the charged lepton Yukawa couplings). This group is broken down to 
U(l) by the Majorana mass term MN2N3. This U(l) group is then only slightly broken by 
the diagonal mass terms ~ AAIm ^ M and by corrections in the Yukawa sector, which can 
be as small as 5inv ~ 0.01 defined in eq. (j3.12p . So, if the existence of slightly broken approx- 
imate symmetry indeed matters, then the inverted hierarchy of neutrino masses with small 
^13 ~ ^inv and small deviation of the angle 8623 ~ Sim from the maximal value is preferred. 

Interestingly, for |e — 1| 

^ *^inv and Sij^Y f the interactions of the heavy neutral lepton 
mass eigenstates with intermediate weak vector bosons are universal and characterized by 
the same mixing angle 

2 _ FV _ matm , . 

Also, both the high temperature baryogenesis and low temperature leptogenesis can take 
place. 

Let us now try to guess what kind of couplings of the dark matter sterile neutrino with 
leptons may lead to some non-trivial symmetries. The phenomenology of DM sterile neutrino 
requires its mass be much smaller than the mass of the singlet fermions responsible for 
baryon asymmetry and that its Yukawa constants are much smaller than those for heavier 
neutral leptons. This leads to a conjecture that the singlet fermion Majorana masses could 
be proportional to their Yukawa couplings, satisfied already for A^2,3 in the construction 
presented above. If true, then the mixing angles of all three sterile leptons with neutrinos 
are the same, and the interaction of them with W and Z bosons exhibits the global SU(3) 
symmetry, which exists for charged leptons. If this hypothesis happens to be correct, the 
mixing angle of DM sterile neutrino is predicted to be 



^Im = ^"'^2''"' = ^1/ = ^ - 2.5 X 10-11 , (8.2) 



corresponding to M ~ 2 GeV, a preferred value leading to the required hierarchy between 
baryon asymmetry and low temperature lepton asymmetry. For this value of the mixing angle 
the mass of DM sterile neutrino is bounded from above by Mi < 8 keV by X-ray observations 
(see the plots presented in ref. [l9]). If the Lyman-a bounds of refs. [44:\ H5] are correct. 
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then Ml > 4 keV (see the discussion in ref. 09]). To exclude or verify this prediction, the 
current X-ray constraints must be improved by a factor of 10. 

Yet another fine-tuning which is necessary for creation of large low-temperature lepton 
asymmetry, is eq. (I2.14p . leading to Scenario Ila for singlet fermion mass difference. Though 
the origin of different parameters even in the Standard Model remains a mystery, it is tempting 
to speculate how this relation, equivalent to 

2{h^h)23v'^ + M{AM^2 + AM33) ~ , (8.3) 



may come from some more fundamental theory. In the i^MSM described by Lagrangian (j2.ip 
the first term in this condition is due to the Higgs condensate while the second is due to 
Majorana masses of singlet fermions and, therefore, they have completely different nature. 
Clearly, a correlation between two independent dimensionfull parameters would be a miracle 
if the z^MSM were the final fundamental theory. This is not so if the mass parameters in the 
z/MSM have the common source, as in the model of [6], where the Higgs boson and the neutral 
fermion masses come from the vacuum expectation value of the nearly conformally coupled 
scalar field singlet with respect to the SM gauge group. In this case the relation (j8.3p turns 
into a connection between the Yukawa coupling constants in the sterile neutrino sector of the 
z/MSM. In fact, all phenomenological and cosmological requirements to the parameters of the 
z/MSM with extra scalar field x can be encoded in a simple Lagrangian, kind of Effective 
Theory of Everything (ETOE), containing just few dimensionless parameters, their powers, 
and one mass scale. It has the form: 

1 X - c 

jCyMSM /^,.MSM[M^o] + 7rr-r(^MX)^ " TT Ni^mijNj + h.c. - V($, x)+Cg, (8.4) 

where the first term is the z/MSM Lagrangian without Higgs potential and with all dimen- 
sionfull parameters (Higgs and Majorana masses) put to zero, the constants /o and /i will 
be specified below. The scalar potential is given by 

V{<f, x) = A (c^tcj. -^^y + I3{x^ - v^f , (8.5) 

where A ~ /? ~ 1/10 are the Higgs and x self-couplings correspondingly, v is the Higgs vev. 
The gravity part is 

,,._(^,.,>.,.)|, (s.) 

where R is the scalar curvature. This is a Lagrangian of "induced gravity" going back to 
refs. [72l [73] (see also ^ in the z^MSM context). The Yukawa couplings hai in eq. (|2.ip are 
written as 

hal = fofajmji , (8.7) 
where faj is an arbitrary complex matrix with elements faj ~ 1 and 



mij = fi 



1 / 

-^0 - [IUMq + transposed 
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with 



\ / 1 

1 + /o a 
10/ V -a 



) 



(8.9) 



where a ~ 1 is a real number. The second term in (j8.8|) is chosen in such a way that eq. (|8.3|) 
is automaticaUy satisfied for any choice of /o, /i and faj- 

The parameter /o = {v/Mp)^ ~ 4 x lO'^, where Mp = (87rGAr)"5 = 2.4 x 10^^ GeV is 
the reduced Planck scale, appears in several places in Lagrangian (j8.4|) . It is fixed from the 
requirement to provide the known Newton constant and the correct phenomenology of singlet 
fermions, as we describe in what follows. 

The parameter /i ~ 6inv ~ 10~^ sets the mass of the singlet fermions (giving baryon and 
lepton asymmetries) in the GeV region, M2_3 ~ fiv. The mass of dark matter neutrino 
Ml ~ /oM2,3 is then in 0(10) keV region, masses of active neutrinos in the fraction of eV 
region, m,^ ~ /oM2,3. In addition, /q makes the Yukawa coupling of sterile neutrino to be 
small compared with the Yukawa couplings of by a factor /o, exactly what is needed 
to produce them in the early universe in right amounts to play the role of dark matter. 
Moreover, A/(2/o) ~ 2 x 10^, appearing in the conformal coupling of the Higgs field to Ricci 
scalar R, leads to inflation producing correct amplitude of primordial fluctuations [13]. The 
field X is very light due to its conformal coupling to gravity (m^^- ~ y/^f^v ~ 10~^ eV), but 
practically decouples from the fields of the i/MSM [74^ [75] (see also [13] ) . In contrast with 
[6], where the terms (18. 6p were not introduced, it plays no role in inflation and in production 
of dark matter sterile neutrinos. 

The author has no idea from were the structures discussed above can be coming from but 
is amazed by some numerical coincidences they uncover. 

9. Conclusions 

In this work we scrutinized the mechanism of leptogenesis via oscillations of light singlet 
fermions and determined the parameter space of the i/MSM which can lead to successful 
baryogenesis. The kinetic processes in the model are quite complicated as they are char- 
acterised by a number of different time scales and by fluctuations (deviations from thermal 
equilibrium) of different nature, interacting with each other. 

The flrst sector includes CP- even deviations from thermal equilibrium in the system of 
almost degenerate singlet fermions. These fluctuations give a "source" term for baryogene- 
sis; creation of lepton asymmetry switches off when these deviations are damped away. The 
kinetic evolution of these fluctuations is governed by four different time scales: two equilibra- 
tion rates for A''2,3, the rate of losing of quantum coherence in oscillations of singlet fermions, 
and the rate of oscillations, related to the mass difference between singlet fermions. In the 
paper we estimated all these time scales (Sec. [5]). We found, in particular, the temperature 
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dependence of the oscillation time, essential for Scenarios I and II for singlet fermion mass 
difference. 

The second sector includes CP- odd deviations from thermal equilibrium in the system 
of singlet fermions and active leptons. There are 7 different essential kinetic time scales 
there. The first 4 are similar to those described above, three others govern the damping of 
asymmetries in different active leptonic flavours. 

We established that the oscillations of must be coherent for effective leptogenesis. 
This is only true if both N2 and are out of thermal equilibrium. In other words, lepton 
asymmetry increases in time till one of the singlet fermions, which interacts more strongly 
with the plasma (A^2 in our notations) enters in thermal equilibrium. After this moment 
the coherence in singlet fermion oscillations is lost, and asymmetries in different quantum 
numbers (which we identified) are damped with the rates, which we determined in Sec. [6j 

We found that baryogenesis may occur in a wide range of singlet lepton masses ranging 
from 140 MeV, allowed by experimental and BBN constraints, to the masses exceeding the 
electroweak scale. An essential requirement is a near degeneracy of a pair of the heavy 
neutral leptons. In addition, the parameter e, characterising the breaking of the U(l) leptonic 
symmetry, cannot be smaller than 7 x 10~^. 

We determined explicitly the CP-violating phase which drives baryogenesis in the model 
and demonstrated that it cannot be expressed only in terms of CP-violating phases of the 
active neutrino mixing matrix. Moreover, we found that the baryon asymmetry is non-zero 
in the limit of small ^13. 

We showed, furthermore, that the i^MSM interactions of singlet fermions may produce 
a significant low temperature lepton asymmetry, being consistent with neutrino oscillation 
experiments and leading to the observed baryon asymmetry of the universe. In a companion 
paper 09], we show that this lepton asymmetry can account for all the dark matter in the 
universe. Thus, the i/MSM without introduction of any new physics or fields such as the 
inflaton may happen to be a correct effective field theory all the way up the Planck scale 
|12j explaining a variety of phenomena that the SM fails to deal with. It is intriguing that 
the production of the baryon asymmetry of the universe and of the dark matter is due to 
essentially the same mechanism, making a step towards understanding why the abundances 
of dark and baryonic matters are roughly the same. 

We also found that large lepton asymmetries in singlet fermions N^, which could have 
been generated above the electroweak scale, may not be in conflict with the observed baryon 
asymmetry and can survive till low temperatures in a specific part of the z^MSM parameter 
space. It corresponds to masses above 140 MeV and small e < 5 x 10~^ and also require 
Scenario Ila for the singlet fermion mass difference. Another possibility is to have singlet 
fermion masses near the pion mass and e in the range 5 x 10--^ < e < 0.01. These regions 
can be explored in kaon experiments and in searches for singlet fermion decays pT]. 

The requirement that the z/MSM produces a lepton asymmetry large enough to speed up the 
dark matter production allows to constrain considerably the parameters of the i^MSM. The 
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most non-trivial requirement is (j2.14p . telling that the zero-temperature difference between 
masses of the physical singlet fermions must be much smaller than the active neutrino mass 
differences. For this choice of parameters the baryon asymmetry is generated at temperatures 
close to the sphaleron freeze-out, T ~ 130 — 175 GeV, and a large lepton asymmetry at 
relatively small temperatures, T = T_ ~ 0.1 — 10 GeV, corresponding to the decoupling of 
singlet fermions from the plasma or to their decays. Later the lepton asymmetry is transferred 
to the dark matter population of sterile neutrinos. The asymmetry generation mechanism 
works for all singlet lepton masses admitted by experimental and BBN constraints discussed 
in [11] and for both types of neutrino mass hierarchies; to produce the low temperature 
lepton asymmetry required for resonant dark matter production the parameter e should be 
large enough, e > 2 x 10~^. Moreover, the requirement of having a much smaller baryon 
asymmetry favours large e ~ 1 and singlet fermion masses in the ©(GeV) range. Particles 
with these properties can be searched for at existing accelerators [H], which is however very 
challenging due to the large value of e, leading to a suppression of their production and to a 
decrease of their decay rates. At the same time, the CP-asymmetry in their decays must be 
at least on the level of few %. 

We speculated on the origin of the necessary fine-tunings in the i^MSM and proposed 
a Lagrangian, containing two dimensionless parameters and their powers, which encodes 
different relations required for the phenomenological success of the model. We found, in 
particular, that the theory with e = 1 and inverted hierarchy of neutrino masses exhibits a 
SU(2) flavour symmetry in the singlet fermion Yukawa sector, broken to U(l) by the Majorana 
mass term. The magnitude of the breaking of this U(l) group is small and is of the order 
4Am^°' — 8 ^ 10~ . In a search of a "maximally symmetric" version of the i/MSM we found 
that it is phenomenologically acceptable to think that the strength of the weak interactions 
of all types of singlet fermions is universal. This conjecture leads to a specific prediction for 
the mixing angle of dark matter sterile neutrino A'^i, potentially testable with the help of 
existing X-ray satellites. 

Finally, a word of warning. All the constraints discussed above are applicable only in 
the case when at temperatures well above the electroweak scale concentrations of all singlet 
leptons are zero. In particular, if the dark matter sterile neutrinos are generated above the 
electroweak scale in right amounts, no generation of large lepton asymmetry is needed below 
the electroweak scale. 
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Appendix A. The rates of singlet fermions production 



We specify in this Appendix the ingredients that went into producing Fig. [3l 

The basic formahsm we follow is that of refs. [HIS]. More precisely, the quantity Y in Fig. [3] 
is given by Eq. (4.8) of ref. [9], while Ym contains the additional weight (go — q)/{qo + q), cf. 
Eq. Km . 

The main difference with respect to the analysis of ref. [S] is that we now consider the 
heavy sterile neutrinos, and that the temperatures are correspondingly higher. This implies 
that the exponentially suppressed 1-loop corrections (Eq. (3.1) of ref. [9]) start to dominate 
over the 2-loop terms (Sec. 3.2 of ref. [9]). More precisely, the main changes to the numerical 
code are as follows: 

• Because of the higher temperatures, the contribution of the bottom quark has been 
added to the 2-loop processes listed in Table 1 of ref. [9]. 

• Once the temperature increases above 20 GeV or so, the treatment of 2-loop effects 
through the Fermi model is no longer justified. Therefore we smoothly switch off the 
2-loop contributions within the range T = (15... 30) GeV. 

• Concerning the 1-loop effects, the graphs to be considered are given in Fig. [2] of the 
present paper (except that we use here a basis where the Majorana mass matrix is 
flavour-diagonal). In the top graph, the particle in the loop can either be a Higgs or a 
Goldstone. In addition, the self-energy of the active neutrino, appearing in the bottom 
graph, depends on the gauge choice (because the active neutrino is off-shell). 

Now, the simplest gauge choice in this context is that of Feynman. Then the real part 
of the active neutrino self-energy (the function b) can be taken directly from ref. [68] 
and the imaginary part from Eq. (3.1) of ref. At the same time, the top graph of 
Fig. [2] amounts to 

^^^^^ ^) = ^ ^^Tr [Q aL Im ^ Higgs an] , (A.l) 

where Im^ Higgs has exactly the form in Eq. (3.1) of ref. [9], with three channels char- 
acterized by PC = l;m-c = mHirnic = ^ua'-, PC = '^'-.mc = mz^mi^ = m^^; and 
PC = 2-mc = mw;'mic = itt-Ic- 

Another possible choice is the unitary gauge. Then the Goldstone contributions can be 
dropped from the top graph, but the active neutrino self-energy needs to be modified. 
We have checked that after the appropriate changes, the numerical results in the two 
gauges differ by an amount which is insignificant on our resolution. 

• Once the temperature increases to several tens of GeV, the evolution of the Higgs 
vacuum expectation value needs to be taken into account. We do this by scaling 
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V2v{T) = 246 GeV^l - r2/r(f, where Tq is fixed through the knowledge that the 
sphaleron freeze-out temperature Tew, where we start our evolution, is characterized 
by \/^v{Tew) — Tew- We choose rriH — 200 GeV and then, according to ref. [62] , 
Tew ~ 175 GeV. All physical particle masses are rescaled by v{T)/v{0). 

Apart from these changes, the numerical techniques used are identical to those in ref. [S]. 
Appendix B. Lower bounds on Yukawa couplings 

In this Appendix we present a lower bound on the following combinations of Yukawa couplings 
which will appear in the analysis of equilibration in the early universe, 

\faa\^={\ha2\^ + \ha3\^) ■ (B.l) 

With the use of (j2.6p one can see that the minimal value of |/q.q.P is simply \[M^]aa\M/v'^. 
The smallest Yukawa couplings correspond to the smallest value of the Majorana neutrino 
mass, which we take to be AI ~ ~ 140 MeV (the mass of the pion is introduced as a useful 
parametrisation) (smaller values would be in conflict with predictions of BBN [69^ [70] and 
experiments devoted to the search of singlet fermions |7H[7l[TT]). Inserting the central values 
for neutrino masses and mixing angles from [56]: Am^^; = 8.0 x 10~^ eV^, Am^^^ = 2.5 x 10~^ 
eV^, 023 = 7r/4, tan^(0i2) = 0.45, ^13 = 0, and choosing the unknown CP-violating phases in 
a way to minimize the Yukawa couplings, we get for the normal hierarchy: 

l/eeP > 1.3 X 10-1^ , > 10-l^ \frr\^ > 10''^ (B.2) 

and for the inverted hierarchy 

l/eeP > 8.8 X 10-1^ , If^^f > 4.4 X 10-l^ \frr\^ > 4.4 X 10'^' . (B.3) 

These numbers change somewhat if the neutrino mixing parameters are varied in the experi- 
mentally admitted ranges. To get a minimal possible value of, say, |/eeP one should take the 
maximal possible atmospheric mass difference (2.7 x 10~^ eV), minimal solar mass difference 
(7.7 X 10~^ eV), minimal 612 — 0.56 and maximal ^13 ~ 0.11, leading to 

l/eeP > 8.4 X 10-18 . (B.4) 

If M > rriTT then the lower bounds are stronger by a factor M/m-,^. 
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